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The  major  purpose  of  this  project  was  to  examine  the  ultra-violet  to 
visible  emission  spectra  produced  by  intense,  relativistic  long  pulse 
electron  beam  interactions  with  gas  and  plasma  environments. 

Two  electron  beam  accelerators  were  used  in  the  interaction 
experiments.  The  first  was  MELBA  (Michigan  Electron  Long  Beam 
Accelerator)  which  operates  at  the  following  parameters: 

Electron  beam  Pulselength:  1  to  4  microseconds 
Beam  Voltage:  -0.7  to  -1 .0  MV 
Beam  Current:  1  -  30  kA,  and 
Energy  per  Pulse:  10-20kJ 

A  second  electron  generator  (Febetron)  with  long-pulse  modules  was 
also  employed: 

Beam  pulselength:  0.4  microseconds  fullwidth, 

Beam  voltage:  -  0.3  to  -0.5  MV 
Beam  current:  1 .4  kA 
Energy  per  pulse:  ~135  Joules 

State  of  the  art  optical  diagnostics  were  used  to  examine  the 
ultra-violet  to  visible  radiation  produced  by  the  interaction  of  the 
electron  beam  with  gases  and  plasmas.  The  light  produced  was  imaged 
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1  onto  the  entrance  slit  (by  means  of  turning  mirrors  and  lenses)  of  either 
1)  a  low  resolution  (0.275  meter)  spectrograph  or  a  high  resolution  (1 
meter)  spectrograph.  The  light  signal  was  detected  with  either  an 
intensified  optical  multichannel  analyzer  in  gates  ranging  between  50ns 
and  4  ps  or  a  photomultiplier  tube.  The  emission  spectroscopy  studies 
concentrated  on  identification  of  the  spectral  lines  from  180  to  800  nm. 
Various  plasma  parameters  were  determined  from  the  optical  emission. 

The  two  types  of  gases  studied  were  rare  gases  (helium,  neon,  argon, 
krypton  and  xenon)  and  air  mixtures  (either  nitrogen  or  dry  bottled  air). 

The  results  of  the  various  studies  are  summarized  below. 
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Two  studies  of  electron  beam  interactions  with  air  were  made.  The 
first  used  MELBA  to  investigate  the  interaction  of  intense  electron  beams 
with  low  pressure,  40  to  100  mTorr,  air.  Light  was  gathered  from  a  point 
45  cm  in  front  of  the  foil  at  an  angle  of  45°.  A  Faraday  cup  was  also  used 
downstream  of  the  beam  to  detect  propagation.  Even  though  the  area  of 
the  beam  sampled  was  quite  small,  the  signal  to  noise  was  adequate 
illustrating  the  fact  that  optical  sensing  can  be  quite  sensitive. 

The  two  major  factors  in  determining  the  amount  of  light  observed 
were  the  pressure  of  the  gas  and  the  amount  of  e-beam  current 
propagation.  Basically,  both  increased  pressure  and  increased  propagatic 
resulted  in  good  signal  to  noise  ratio.  However  propagation  generally 
occurs  for  a  pressure  window  of  50-500  mTorr  and  instabilities  result  as 
the  pressure  increases  above  100  mTorr.  Basically  the  best  optical 
signatures  occurred  for  pressures  between  40  and  75  mTorr.  Note  that  foi 
all  pressures  the  air  emission  spectra  looked  similar;  The  spectra  of  dry 
air  showed  strong  emission  from  the  first  negative  band  of  N2+ 
(B2E->X2Z).  Atomic  neutral  and  ion  lines  were  never  observed,  nor  were 
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t  impurities  such  as  hydrogen  or  oxygen. 

The  second  set  of  studies  were  performed  using  the  Febetron.  The 
pressures  of  these  experiments  ranged  from  about  1  Torr  to  100  Torr.  The 
spectra  looked  similar  to  the  low  pressure  cases  observed  with  MELBA  but 
strong  emission  from  the  second  positive  band  of  neutral  N2  (C3n->B3n) 
were  also  observed.  Again  no  other  components  of  air  (including  atomic 
nitrogen)  were  detected. 

The  vibrational  temperatures  of  the  observed  vibrational  bands  were 
determined  from  molecular  Boltzmann  plots  to  be  about  0.33  eV  for  the  N2 
emission  and  0.15  eV  for  the  N2+  emission.  These  temperatures  fall 
between  the  electron  temperature  of  1  -  2  eV  estimated  from  the  rare  gas 
atomic  Boltzmann  plots  (discussed  below)  and  the  gas  temperature  of 
0.025  eV  estimated  from  heating  and  thermal  equilibration  calculations. 

The  vibrational  temperature  is  not  necessarily  equal  to  the  electron 
temperature  because  vibrational  populations  require  more  collisions  to 
equilibrate  than  do  free  electron  populations.  The  vibrational  temperature 
is  smaller  than  the  plasma  electron  temperature  because  the  massive 
nitrogen  atoms  resist  strong  vibrations,  but  is  larger  than  the  gas 
temperature  because  the  plasma  does  transfer  energy  to  the  molecule's 
vibrational  mode. 

Note  that  studies  were  also  made  of  purely  microwave  excited  air 
discharges  in  an  attempt  to  distinguish  the  optical  emission  due  to  the 
electron  beam  alone  and  the  microwaves  generated  by  the  beam.  The 
optical  emission  of  microwave  (2.45  GHz)  discharges  was  similar  to  the 
electron  beam  produced  discharges  in  that  emission  from  N2  was  observed 
with  the  absence  of  atomic  lines.  In  the  case  of  the  pure  microwave 
discharges  however,  the  N2+  lines  were  much  lower  in  intensity  than  the 
electron  beam  produced  discharges  leading  to  the  conclusion  that  the 
electron  density  of  the  e-beam  discharges  was  much  higher.  In  addition,  a 
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t  long  lived  continuum  (~  1  second)  was  observed  in  the  microwave 
discharge  due  to  NO  recombination  with  O  atoms.  The  vibrational 
temperatures  obtained  were  slightly  higher  than  the  e-beam  plasma  alone. 
They  ranged  from  0.3  to  0.7  eV  for  the  Ng  emission.  Basically  the  major 
difference  in  distinguishing  the  electron  beam  produced  optical  emission 
and  microwave  produced  optical  emission  was  the  NO2  afterglow. 

Electron  Beam-Generated  Emissions  from  Rare  Gases 

A  comprehensive  study  of  optical  emission  produced  by  electron 
beams  from  the  Febetron  was  made.  There  were  several  characteristics 
which  were  common  to  all  of  the  gases  studied  with  the  Febetron:  1)  In 
all  cases,  the  light  emission  fell  to  zero  when  the  beam  pulse  ended,  in 
contrast  to  many  REB-pumped  plasma  experiments  where  the  fluorescence 
of  the  decaying  plasma  is  measurable  for  tens  of  ns  to  several  jis  after  the 
beam  ends.  This  indicates  that  the  beam  electrons  are  very  important  in 
the  population  of  excited  states  in  this  experiment.  2)  The  time  histories 
of  the  plasma  showed  little  variation  as  the  plasma  evolved.  The  relative 
intensities  of  the  emitted  lines  were  roughly  constant  over  time, 
indicating  that  the  plasma  reached  steady  state  within  the  first  50  ns  of 
the  beam  pulse.  This  is  is  consistent  with  heating  arguements  which 
indicate  that  both  the  beam  and  return  current  heating  reach  steady  state 
after  the  50  ns  risetime  of  the  beam  current.  3)  Rare  gas  excimer 
emission  was  not  detected  due  to  the  subatmospheric  pressures  used.  The 
VUV  excimer  laser  lines  were  not  examined,  but  several  less  prominent 
visible  and  ultraviolet  lines  would  have  been  detected  if  excimer  emission 
was  significant.  This  is  in  agreement  with  the  chemical  kinetics  model 
which  predicts  excimer  emission  to  be  negligible.  Also,  unusual  molecular 
emission  due  to  molecular  impurities  was  not  observed.  4)  All  of  the  rare 
gases  exhibited  beam-dominated  excitation  regimes  at  low  gas  pressures 
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t  '  where  the  intensity  of  the  emitted  light  was  linear  with  gas  pressure. 

The  intensities  fell  off  at  pressures  above  400  Torr,  possibly  due  to 
self-absorption  of  the  emitted  photons. 

The  spectra  emitted  by  the  rare  gases  followed  a  trend  that  scaled 
with  ionization  potential  in  that  singly-ionized  lines  dominated  the 
spectra  of  the  more  weakly-binding  gases  (Ar,  Kr,  and  Xe),  while  neutral 
lines  dominated  the  spectra  of  the  more  strongly-binding  gases  (He  and 
Ne).  Doubly-ionized  lines  were  most  prevalent  in  xenon,  which  has  the 
lowest  binding  energy  of  the  rare  gases,  and  were  not  observed  in  gases 
lighter  than  krypton.  Rare  gas  excimer  emission  was  not  observed  due  to 
the  subatmospheric  pressures  used  which  inhibited  the  collisional 
formation  of  rare  gas  molecules.  The  pressure  dependence  of  the  rare  gas 
emission  shows  intensity  to  be  roughly  linear  with  gas  pressure  between 
25  to  300  Torr.  The  light  intensity  falls  off  about  300-700  Torr,  possibly 
due  to  self-absorption  in  the  relatively  denser  plasmas. 

The  optical  emission  was  modeled  by  extending  the  Local 
Thermodynamic  Equilibrium  (LTE)  model  to  account  for  the  presence  of 
fast  beam  electrons.  These  energetic  electrons  can  populate  high-lying 
energy  states  far  above  their  expected  LTE  levels,  resulting  in  intense 
light  emission  from  these  states  as  they  spontaneously  decay.  Indeed,  the 
beam  electrons  are  almost  entirely  responsible  for  the  production  of 
excited  ionic  states  in  this  experiment  because  the  plasma  densities  of 
1015-1016  cm'3  calculated  for  this  experiment  by  Brake  and  Repetti  are 
too  high  to  populate  the  rare  gas  ionic  states  through  normal  LTE 
collisional  processes.  Assuming  that  the  ionic  states  are  populated  via 
excitation  from  the  neutral  ground  state  in  this  experiment,  it  is  possible 
to  correct  for  the  beam  enhancement  and  put  both  ionic  and  neutral  lines 
on  an  atomic  Boltzmann  plot  to  retrieve  the  excitation  temperature  of  the 
plasma. 
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<  The  excitation  temperature  equals  the  plasma  electron  temperature 

if  the  plasma  has  achieved  complete  local  thermodynamic  equilibrium. 

The  plasma  densities  in  this  experiment  are  not  large  enough  to  achieve 
full  LTE,  but  are  large  enough  to  achieve  partial  LTE  where  the  populations 
of  the  excited  states  observed  in  emission  are  close  to  their  LTE  values. 

The  excitation  temperatures  obtained  by  the  modified  Boltzmann  plots  fell 
in  the  range  of  1 .5  ±  0.7  eV  for  the  rare  gases,  increasing  as  atomic 
weight  decreased. 

The  most  probable  reason  that  the  effect  of  e-beam  enhancement  has 
not  been  reported  previously  is  the  unique  beam-plasma  regime  studied  in 
this  experiment.  The  Febetron  electron  beam  was  not  intense  and  sprayed 
out  from  the  anode  foil  into  a  40°  half-angle  cone  as  int  entered  the 
plasma.  The  beam-to-plasma  density  ratio  of  10-4  to  10‘6  was  very  low, 
reducing  the  return  current  and  instability  heating  of  the  plasma.  This 
kept  the  plasma  secondary  electron  temperature  relatively  low  at 
approximately  1  eV,  which  was  not  enough  to  appreciably  populate  the 
observed  ionic  states  at  the  existing  plasma  densities.  Under  these 
conditions,  the  beam  excitation  effect  was  clearly  evident  in  the  emission 
spectra. 

Other  factors  contributed  to  the  unique  nature  of  this  expeirment. 

The  300  ns  duration  of  the  electron  beam  was  not  long  enough  for  the 
electrons  to  appreciably  heat  the  ions,  but  was  long  enough  to  achieve  a 
steady  state  equilibrium  in  the  plasma.  The  subatmospheric  pressures 
were  low  enough  to  discourage  excimer  formation  but  were  high  enough  to 
achieve  partial  Local  Thermodynamic  Equilibrium  for  the  observed 
transitions  and  also  kept  the  secondary  plasma  electron  temperature  low 
through  increased  collisions.  The  plasmas  were  optically  thin  at  low  to 
intermediate  pressure  which  precluded  consideration  of  radiative 
excitation  and  photon  transport.  Finally,  the  optical  multichannel 
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anlayzer  allowed  examination  of  a  wide  range  of  wavelengths.  All  of 
these  factors  help  explain  the  differences  between  the  optical  emission 
observed  herein  and  the  results  of  other  investigations. 

Studies  of  helium  made  with  MELBA  at  pressure  between  75  and  227 
mTorr  showed  emission  due  only  to  neutral  helium.  Electron  temperatures 
estimated  from  electronic  temperatures  determined  from  ratios  of  line 
intensities  yielded  a  temperature  of  about  0.5  eV. 
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Abstract — Experiments  have  been  performed  demonstrating  the  ef¬ 
ficient  transj— • ;  of  long-pulse,  high-current  electron  beams  in  well- 
diagnosed,  preformed  plasma  channels.  Channels  were  created  by  the 
low-energy  electron-beam  ionization  of  low-pressure  monatomic  gases. 
These  transport  efficiencies  are  comparable  to  the  optimum  values  ob¬ 
tained  for  similar  putsetengtbs  in  previous  long-pulse,  electron-beam 
transport  experiments  in  laser-pre ionized  diethylaniline  gas.  Experi¬ 
mental  results  show  the  development  of  both  transverse  current  cen¬ 
troid  oscillations  characteristic  of  the  ion  hose  instability  and  the  mi¬ 
crowave  emission  indicative  of  relativistic  electron  beam-plasma 
electron  two-stream  instabiliiy.  Significant  microwave  emission  indi¬ 
cating  rapid  growth  of  the  two-stream  instability  has  been  observed  in 
the  5-band  when  the  space-charge  neutralization  fraction  /,  exceeds 
unity. 

I.  Introduction 

PULSED,  high-power  relativistic  electron-beam  (REB) 
technology  has  significant  applications  to  advanced 
accelerators,  the  generation  of  high-power  microwaves 
[1],  laser  excitation  [2],  and  free  electron  laser  develop¬ 
ment  [3].  In  particular,  recent  advances  in  pulsed  power 
technology  [4]  have  made  available  electron-beam  gen¬ 
erators  capable  of  producing  high-current  REB's  of  long 
duration  (  >  l  ns).  The  benefits  to  these  technologies  of 
large  delivered  energies  ( >  10  kJ)  using  long-pulse  elec¬ 
tron  beams  may  be  significant. 

Important  to  these  applications  is  the  stable  transport  of 
the  electron  beam  over  distances  on  the  order  of  meters. 
Stable  propagation  of  high-current  REB’s  in  the  ion  focus 
regime  (IFR)  is  based  on  the  generation  of  a  well-defined 
plasma  background,  created  by  REB-neutral  gas  colli- 
sional  ionization  or  by  preionizing  a  channel  in  low-pres¬ 
sure  gas;  plasma  electrons  are  immediately  expelled  from 
the  channel  by  the  space  charge  of  the  electron  beam.  The 
field  of  the  remaining  plasma  ions  reduces  the  space- 
charge  potential  of  the  electron  beam  enough  to  allow  the 
electrostatic  repuli  ive  and  self-magnetic  field  forces  of  the 
beam  to  cancel,  allowing  stable  propagation. 
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We  present  results  of  experiments  to  investigate  the 
physics  involved  in  the  interaction  between  a  long-pulse 
high-current  REB  with  preformed  IFR  plasma  channels. 
The  behavior  of  the  electron  beam-channel  plasma  inter¬ 
actions  may  not  be  as  straightforward  as  that  for  short 
pulse  electron  beams,  since  the  timescales  may  now  be¬ 
come  long  enough  for  the  dynamics  of  the  channel  ions 
to  become  significant. 

Short-pulse  electron-beam  transport  in  preformed  IFR 
plasma  channels  was  initially  investigated  elsewhere  using 
channels  preformed  in  organic  gases  (benzene)  [S],  [6]  by 
multiphoton  UV  ionization.  Previous  experiments  at  the 
University  of  Michigan  extended  this  technique  to  the  in¬ 
vestigation  of  the  characteristics  of  long-pulse  (  *  1  #zs) 
electron-beam  transport  [7],  [8]  in  channels  formed  using 
a  KrF  excimer  laser  to  preionize  a  channel  in  diethylani¬ 
line  (DEA). 

Photoionization  of  organic  molecules  by  UV  laser  has 
drawbacks,  however  [9],  for  accelerator  transport  sys¬ 
tems.  The  organic  compounds  can  be  difficult  to  handle, 
and  fragmentation  of  the  large  molecules  can  complicate 
the  channel  physics.  Additionally,  to  prevent  the  channels 
from  becoming  overneutralized  due  to  beam-induced  ion¬ 
ization  during  the  beam  pulse,  very  low  gas  pressures  and 
high-power  UV  lasers  are  required  to  ionize  the  gas  to  the 
required  level.  This  is  particularly  significant  for  the 
transport  of  long-pulse  beams. 

The  experiments  reported  here  utilize  low-energy  elec¬ 
tron  impact  ionization  to  overcome  the  limitations  of  long- 
pulse  laser  guiding.  Since  electron-impact  ionization  cross 
sections  are  large,  the  background  gas  is  not  restricted  to 
the  use  of  organic  molecules,  and  monatomic  background 
gases  may  be  used  which  simplify  the  ion  channel  physics 
[  10]— [  12].  Typically,  a  weak  longitudinal  magnetic  field 
(  <  100  G)  is  applied  to  confine  the  low-energy  electron 
beam  and  produce  well-defined  plasma  channels.  Since 
monatomic  gases  have  a  relatively  long  recombination 
lifetime,  the  channel  ionization  can  build  up  over  a  long 
period  of  time  (several  hundred  microseconds)  with  a  very 
low-current,  low-energy  beam  [10]. 

II.  Plasma  Channel  Formation 

These  experiments  investigate  long-pulse  (350  ns)  REB 
propagation  in  pure-mass  monatomic  ion  channels. 
Plasma  channels  have  been  formed  by  low-energy  elec- 
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current  centroid 
displacement  coils  and 
net  current  monitor 


Fig.  1.  Experimental  configuration  for  long-pulse  REB  injection  into 
pieionized  plasma  channels. 


Fig.  2.  Channel  plasma  density  on  axis  54  cm  from  the  filament  versus 
background  gas  pressure.  Discharge  parameters:  800  mA,  -300  V.  and 
confining  axial  magnetic  field  of  100  G.  (A)  —  neon:  (o)  ■»  argon:  (°) 
«=  krypton. 

tron  impact  ionization  of  noble  gases  (He.  Ne,  Ar,  Kr, 
and  Xe).  The  experimental  configuration  is  shown  in  Fig. 

1 .  The  plasma  source  is  a  heated  tungsten  filament  on  axis, 
dc-biased  to  —300  V  with  respect  to  the  walls  of  the  trans¬ 
port  tube  (ground).  The  filament  is  immersed  in  the  mirror 
region  of  a  weak  (SO- 100  G)  pulsed  magnetic  field  (7-ms 
risetime).  One  end  of  the  mirror  is  matched  to  the  sole- 
noidal  field  of  the  transport  tube  ( 1-m  long,  15-cm  diam ). 
The  opposite  end  of  the  mirror  has  a  larger  magnetic  mir¬ 
ror  ratio 'to  preferentially  allow  source-plasma  electrons 
to  be  directed  into  the  transport  tube.  Axial  and  radial 
plasma  density  profiles  have  been  measured,  as  well  as 
the  temporal  evolution  of  these  channels  using  Langmuir 
probe  techniques.  The  channel  plasma  density  is  adjust¬ 
able  up  to  about  4  x  10"  cm-3  by  varying  the  back¬ 
ground  gas  pressure,  the  source  plasma  discharge  current,, 
or  the  confining  axial  magnetic  field.  The  maximum  chan¬ 
nel  plasma  densities  obtained  in  argon,  neon,  and  krypton 
are  shown  as  a  function  of  gas  pressure  in  Fig.  2.  Mea¬ 
surements  of  the  radial  channel  Kiasma-density  profiles 
indicate  that  they  are  approximately  Gaussian  with  a 
FWHM  radius  from  1.0  to  l.S  cm.  Axial  plasma  density 
measurements  indicate  uniformity  to  within  20%  down  the 
length  of  the  transport  tube  for  the  gas  pressures  utilized 
in  the  experiments. 


The  channel  plasma  density  is  carefully  monitored  with 
a  Langmuir  probe  by  pulsing  the  channel  separately  be¬ 
fore  injection  of  the  REB.  The  temporal  evolution  of  the 
channel  plasma  density  with  the  confining  axial  magnetic 
field  is  recorded  and  used  to  determine  the  time  at  which 
to  inject  the  electron  beam  on  the  subsequent  pulse.  Since 
the  plasma  channels  have  been  found  to  be  very  repro¬ 
ducible,  this  provides  an  accurate  determination  of  the 
channel  plasma  density  just  prior  to  electron-beam  injec¬ 
tion. 

III.  Electron-Beam  Injection 

Relativistic  electron  beams  are  generated  by  a  long- 
pulse  Febetron1  pulser,  operating  with  the  parameters: 
voltage  =  -250  to  -400  kV;  diode  current  =  1.2  kA; 
and  pylse  length  =  approximately  380  ns.  Electron  beams 
are  generated  from  a  field-emission  cathode-emitting  sur¬ 
face  consisting  of  a  2-cm  square  array  of  carbon  fiber  bun¬ 
dles  spaced  4-mm  apart.  This  emitting  surface  is  set  into 
a  Pierce-type  cathode  geometry  with  a  diode  gap  spacing 
of  2  cm.  A  6- jim-thick  aluminized  Mylar  foil  separates 
the  diode  chamber  from  the  transport  region;  aluminiza¬ 
tion  is  always  kept  facing  the  cathode.  This  thin  Mylar 
foil  is  used  to  minimize  the  scattering  of  the  beam  ex¬ 
tracted  from  the  diode.  After  passing  through  the  anode 
foil  the  electron  beam  is  apertured  by  a  1 -cm -radius  hole 
in  a  (1-mm  thick)  carbon  plate  located  2-cm  downstream 
from  the  anode  foil.  By  keeping  the  aperture  diameter  at 
2  cm,  the  electrons  extracted  from  the  periphery  of  the 
cathode-emitting  surface,  typically  with  large  angular  di¬ 
vergence,  are  sc  raped -off  before  injection  into  the  IFR 
plasma  channels.  Typical  injected  beam  currents  are  ap¬ 
proximately  200  A. 

The  carbon  aperture  plate  can  be  interchanged  with  a 
range-thick  (1  mm),  carbon  pinhole  mask  for  character¬ 
ization  of  the  emittance  of  the  electron  beam  injected  into 
the  transport  tube.  The  pinholes  have  a  1-mm  diam  and 
are  placed  in  a  square  array  6-mm  apart.  A  3.2-mm-thick 
plate-glass  Cerenkov  imaging  screen  is  located  4 -cm 
downstream  from  the  mask.  The  diode  side  of  the  plate- 
glass  Cerenkov  emitter  is  heavily  coated  with  carbon 
aerosol  (Aerodag)  to  block  the  visible  light  emission  from 
the  diode  region.  The  time-integrated  Cerenkov  light  from 
the  beam  electrons  is  viewed  through  a  Lucite  window 
and  a  plane  mirror  set  at  about  45°  to  the  electron-beam 
axis.  Previous  experimental  evaluation  of  Cerenkov  con¬ 
verters  for  measuring  current  density  distributions  in  high- 
current  REB  experiments  have  shown  commercial  plate 
glass  to  have  excellent  time  response  and  correlation  of 
light  output  to  electron  current  density  [13]. 

Time-resolved  transverse  motion  of  the  electron  beam 
is  monitored  with  an  array  of  fi*-dot  (Mimov)  coils  [14], 
[15]  located  around  the  circumference  of  the  transport 
tube.  Location  of  the  current  centroid  (in  x-y  space)  is 
determined  by  two  pairs  of  differencing  coils  (at  right  an¬ 
gles  to  each  other),  together  with  a  measurement  of  the 
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net  current  at  the  probe  location  using  a  Rogowski  coil. 
Two  such  probes  are  used  in  these  experiments.  One 
probe,  located  3-cm  downstream  from  the  carbon  aperture 
plate,  provides  characterization  of  the  electron  beam  at 
injection.  The  second  probe  is  located  95-cm  downstream 
from  the  injection  point  and  monitors  the  behavior  of  the 
transported  beam. 

The  transported  beam  current  is  measured  by  a  fast  rise¬ 
time  Faraday  cup  f  16]  comprised  of  a  vacuum-tight  as¬ 
sembly,  which  is  separately  evacuated  by  a  dedicated  me¬ 
chanical  pump.  It  measures  the  transported  beam  current 
within  2.54  cm  of  the  axis  of  the  transport  tube.  A  25.4- 
pm  thick  titanium  foil  discriminates  against  secondary 
(plasma)  electrons,  filtering  out  electrons  below  70-keV 
energy.  The  transported  electron-beam  radial  current  pro¬ 
file  was  inferred,  with  the  Faraday  cup  removed,  from  the 
damage  pattern  observed  on  a  Lucite  window  at  the  end 
of  the  transport  tube. 

Radio  frequency  (RF)  and  microwave  radiation  spectra 
are  also  monitored  to  investigate  the  initiation  of  stream¬ 
ing  instabilities.  Input  signals  are  derived  from  a  broad¬ 
band  (single  turn)  Ba-dot  loop.  RF  signals  are  fed  to  a 
spectrum  analyzer,  consisting  of  an  array  of  bandpass  fil¬ 
ters  and  diode  detectors,  which  separates  the  signal  into 
four  frequency  bands:  45  to  105  MHz;  138  to  184  MHz; 
210  to  265  MHz;  and  313  to  373  MHz.  Microwave  emis¬ 
sions  are  monitored  using  a  series  of  successively  higher 
frequency  directional  couplers  and  diode  detectors.  The 
signal  is  separated  into:  5-band  (2. 1-4.3  GHz);  7-band 
(4.4-6.6  GHz);  X-band  (6.6-14.2  GHz);  A-band  (14.3- 
30  GHz);  and  Ad-band  (2 1 . 1  -40  GHz). 

The  diode  region  is  pumped  down  to  about  2  x  10-5  T 
by  an  oil  diffusion  pump.  The  transport  tube  is  separately 
evacuated  to  a  base  pressure  of  about  4  x  10-6  T  before 
backfilling  with  the  desired  channel  gas.  Due  to  the  low 
background  gas  pressures  required,  the  diffusion  pump  on 
the  transport  tube  pumps  continuously.  This  dynamic-fill 
technique  has  proven  to  be  quite  reliable  and  very  stable. 

IV.  Experimental  Results  and  Discussion 

Electron-beam  transport  results  have  been  obtained 
using  pre-ionized,  low-density  plasma  channels  pre¬ 
formed  in  low-pressure  neon,  argon,  krypton,  and  xenon 
gas.  For  the  gas  pressures  utilized  in  these  experiments, 
relativistic  electron  impact  ionization  (An,)  is  insignifi¬ 
cant,  since  An,/n *»  0.01  after  350  ns.  Beam  injection 
into  the  transport  tube  without  energizing  the  tungsten  fil¬ 
ament  resulted  in  virtually  no  signal  detected  from  the  net 
current  monitor  located  at  the  end  of  the  1-m-long  trans¬ 
port  tube.  This  provides  confirmation  that  (a)  The  back¬ 
ground  gas  pressures  are  low  enough  that  REB-neutral  gas 
collisional  ionization  does  not  play  a  significant  role  in 
channel  formation;  and  (b)  the  low  axial  magnetic  field 
used  to  confine  the  channel  plasma  has  very  little  influ¬ 
ence  on  the  beam  dynamics. 

The  REB  at  injection  is  characterized  by  a  beam  current 
of  approximately  200  A  with  a  pulse  length  of  380  ns. 
Time-integrated  emittance  measurements  characterizing 
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the  beam  at  injection  indicate  the  initial  angular  diver¬ 
gence  of  the  beam  to  be  approximately  0.35  rad  ( *  20°) 
on  the  envelope.  The  corresponding  unnormalized  emit¬ 
tance  is  determined  to  be  approximately  200  cm  -  mrad. 
Since  the  Cerenkov  threshold  energy  in  plate  glass  is  ap¬ 
proximately  130  keV,  these  time-integrated  emittance  es¬ 
timates  represent  the  first  325  ns  of  the  electron  beam 
pulse. 


A.  Transport  Efficiency 

Efficient  IFR  electron-beam  transport  is  expected  tc  oc¬ 
cur  for  space-charge  neutralization  fraction  values  in  the 
range  fe  &  1  /y2  (where  fe  is  the  ratio  of  the  channel  ion 
density  to  the  electron-beam  density,  fr  -  n,//irt).  For  a 
cold  (i.e.,  zero  emittance)  beam  in  radial  force  equilib¬ 
rium,  the  IFR  condition  is  fe  =  l/y2,  which  for  these 
experiments  (y  *  1.6)  occurs  for/,  =  0.4.  In  practice, 
however,  values  of  fe  greater  than  1  /y2  are  required  to 
compensate  for  nonzero  beam-emittance  effects.  Fig.  3 
shows  a  comparison  between  the  transported  beam  current 
measured  by  the  Faraday  cup  (2. 54 -cm  radius)  and  the 
transported  net  current  measured  by  a  Rogowski  coil  (7.5- 
cm  radius)  for  the  channel  providing  a  space-charge  neu¬ 
tralization  fraction  fe  of  0.5.  This  figure  indicates  that  not 
all  of  the  beam  current  transported  is  collected  by  the  Far¬ 
aday  cup.  This  result  and  the  corresponding  limitations  to 
transport  efficiency  may  be  interpreted  in  terms  of  the 
nonzero  emittance  of  the  REB. 

The  equilibrium  beam  radius  of  the  REB  in  the  IFR 
plasma  channel  may  be  determined  from  a  simple  enve¬ 
lope  equation  analysis  [17]  as: 


*  <2  4.  f?  7 

beq  -2/„(l  -fr-  02) 


(1) 


where  e  is  the  unnormalized  beam  emittance,  70  = 
4ir€omc3/e,  0  =  v./c,  and  y  =  (I  -  02)~'^2.  For  the 
measured,  time-integrated,  injected  electron-beam  emit¬ 
tance  of  2.0  x  10'3  m  •  rad  and  injected  currents  of  200 
A,  the  equilibrium  beam  radius  ranges  from  3.2  cm  for  fe 
—  0.5  to  1.5  cm  for  fc  —  1.0.  For  f,  values  of  0.5,  the 
injected  beam  relaxes  to  an  equilibrium  radius  larger  than 
the  radius  of  the  collecting  surface  of  the  Faraday  cup.  As 
a  result,  not  all  of  the  beam  current  is  collected,  consistent 
with  the  observations.  Radial  profile  measurements  in¬ 
ferred  from  the  damage  pattern  obtained  on  a  Lucite  wit¬ 
ness  plate  produced  with  a  space-charge  neutralization 
fraction  in  argon  of  0.9  are  in  good  agreement  with  the 
radius  predicted  by  (1).  Also  consistent  with  this  model 
is  the  observation  that  as  f,  approaches  unity,  the  trans¬ 
ported  beam  signal  from  the  Faraday  cup  and  the  net  cur¬ 
rent  measurement  are  in  very  good  agreement. 

The  peak  instantaneous  current  transport  efficiency  (de¬ 
fined  as  the  ratio  of  the  transported  current  measured  by 
the  Faraday  cup  to  injected  current)  was  as  high  76%  in 
these  experiments.  For  long-pulse  electron-beam  trans¬ 
port,  an  important  figure  of  merit  characterizing  the  trans¬ 
port  technique  under  investigation  is  the  charge  transport 


MILLER  A*Nl>  GILGENBACH:  TRANSPORT  ME  LONG-PULSE.  HIGH-CURRENT  ELECTRON  BEAMS 


661 


Fig.  3.  Comparison  of  transported  beam  current  (collected  in  the  2.54-cm 
radius  Faraday  cup)  and  transported  net  current  (collected  in  the  entire 
transport  tube  cross  section.  7.5-ctn  radius)  to  the  injected  current  for  /, 
-0.5. 

efficiency.  We  define  the  charge  transport  efficiency  as  the 
ratio  of  the  transported  electron-beam  charge  (as  mea¬ 
sured  by  the  Faraday  cup)  to  the  injected  electron-beam 
chaige.  Summaries  of  the  experimentally  measured  peak 
instantaneous  current  transport  efficiency  and  charge 
transport  efficiency  are  shown  in  Fig.  4(a)  and  (b).  These 
plots  show  the  transport  efficiencies  of  the  preionized 
plasma  channels,  in  all  four  gases,  for  fr  values  in  the 
range  0  to  1.6.  As  expected,  current  (charge)  transport 
efficiencies  greater  than  50%  (40%)  begin  around  f,  = 
0.35,  in  very  good  agreement  with  the  initial  IFR  condi¬ 
tions  of  1  /V  *  0.4.  The  peak  current  transport  efficiency 
appears  to  saturate  at  about  65-70%  (peak  values  of  76%) 
for  ft  values  greater  than  0.9.  Similarly,  the  charge  trans¬ 
port  efficiency  saturates  at  about  50%  (peak  values  of 
56%)  for  ft  values  greater  than  0.75.  This  limit  in  trans¬ 
port  efficiency  is  most  likely  due  to  the  electron  beam  not 
being  well  matched  at  injection  to  the  equilibrium  radius 
of  the  REB  in  the  IFR  channel. 

B.  Beam  Dynamics 

Fig.  5  presents  electron-beam  transport  data  for  beam 
injection  into  an  argon  channel  with  a  space-charge  neu¬ 
tralization  fraction  fr  *»  1.2.  (Note:  ft  values  greater  than 
unity  indicate  the  presence  of  excess  plasma  electrons  in 
the  channel  that  cannot  be  expelled  by  the  space  charge 
of  the  REB.)  Efficient  propagation  occurs  over  the  entire 
pulse  length.  Peak  instantaneous  current  transport  effi¬ 
ciency  is  70%  for  this  pulse,  while  the  total  fraction  of 
electron-beam  charge  transported  within  the  Faraday  cup 
diameter  was  44%.  Plasma  current  effects,  including  the 
return  plasma  current  on  the  injected  net  current  monitor 
and  possible  current  multiplication  on  the  downstream  net 
current  monitor,  were  not  observed  until/,  exceeded  1.2. 

Two  important  features  are  apparent  in  the  data  of  Fig. 
5.  First,  the  downstream  current  centroid  displacement 
probe  indicates  a  sinusoidal  oscillation  in  the  y  (arbitrary) 
plane.  Fig.  6  presents  the  approximate  x-y  motion  un¬ 
folded  from  the  data  of  Fig.  5.  Transverse  oscillation  of 
the  beam  current  centroid  is  clearly  illustrated.  These  re¬ 
sults  are  consistent  with  the  transverse  motion  of  the  beam 
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(b) 

Fig.  4.  (a)  Peak  instantaneous  current  transport  efficiency:  and  (b)  charge 
transport  efficiency  versus /,:  (®)  —  argon  at  5  x  10~*  T:  (°)  -  neon  at 
3  x  10‘4  T;  (A)  -  krypton  at  5  x  10'*  T:  and  (O)  -  xenon  at  3  x 
10'*  T. 


WO  na /*» 

Fig.  5.  Experimental  data  showing  transport  results  for  long -pulse  elec¬ 
tron-beam  injection  into  an  argon  plasma  channel  with  fr  -  1.2.  The 
transported  net  current  monitor  is  low-pass  Altered  at  20  MHz. 
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Fig.  6.  Unfolded  x-y  (arbitrary)  beam  centroid  motion  from  the  data  of 
Fig.  5.  The  data  points  are  spaced  5-ns  apart. 


due  to  growth  of  the  ion  hose  instability  [8],  [12],  [18]- 
[20].  For  these  experimental  conditions  the  theoretical 
frequency  for  ion  hose  motion  [19]  is  about  4.6  MHz, 
which  is  in  reasonable  agreement  with  the  experimentally 
observed  oscillation  frequency  of  2.5  MHz.  Data  taken 
over  a  range  of  fe  from  0  to  1.6  showed  increasing  real 
frequency  and  maximum  excursion  with  increasing  values 
of  fe.  This  is  qualitatively  consistent  with  the  expected 
scaling  [12]  of  the  growth  of  the  ion  hose  instability  (G 
ot  ft  )•  The  development  of  these  oscillations  was  also 
investigated  as  a  function  of  the  channel  ion  mass.  Nearly 
complete  oscillations  were  observed  in  the  neon  {A  =  20 
amu)  and  argon  (/4  =  40  amu)  data;  however,  complete 
oscillations  were  not  observed  in  krypton  (A  *=  84  amu) 
or  xenon  (A  =  132  amu).  This  scaling  with  the  channel 
ion  mass  is  also  qualitatively  as  expected  for  the  ion  hose 
instability  ( G  «  [12]). 

Secondly,  examination  of  the  microwave  emission  sig¬ 
nal  received  by  the  RF  probe  in  the  data  of  Fig.  5  illus¬ 
trates  a  significant  frequency  component  in  the  5-band. 
The  evolution  of  the  microwave  signal  with  the  channel 
plasma  density  is  shown  in  Fig.  7(a)  and  (b).  The  data 
shows  clear  growth  of  the  5-band  signal  with  increasing 
ft,  becoming  significant  at  ft  *  1.1.  This  radiation  ap¬ 
parently  originates  from  the  REB-plasma  electron  two- 
stream  instability  [21].  This  instability  can  occur  when 
the  channel  plasma  density  becomes  great  enough  that  ex¬ 
cess  plasma  electrons  remain  in  the  channel  and  interact 
with  the  beam  electrons.  Typically,  this  streaming  insta¬ 
bility  occurs  for  fr  greater  than  unity  [6].  For  f,  values  on 
the  order  of  1.2,  the  electron  plasma  frequency  is  approx¬ 
imately  1.3  GHz,  near  the  observed  signal  frequency 
range  in  the  5-band.  Microwave  signals  as  shown  in  Figs. 
5  and  7  were  only  observed  on  the  5-band  detector.  No 
such  signal  was  obtained  for  frequencies  above  4.3  GHz 
nor  below  343  MHz.  No  signal  was  observed  in  the 
5-band  frequency  range  for  f,  less. than  0.8  or  0.9  in  any 
gas. 

V.  Summary 

The  data  demonstrates  that  the  large  value  of  emittance 
characterizing  the  REB  at  injection  requires  that  space- 
charge  neutralization  fractions  f,  close  to  (but  less  than) 
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Fig.  7.  (a)  5-band  (2. 1-4.3  GHz)  microwave  output  power  detected  for 
plasma  channels  preformed  in  argon;  injected  beam-current  pulse  is 
shown  for  comparison,  (b)  Plot  of  peak  5-band  microwave  output  power 
versus  (o)  *  argon  at  5  x  10"*  T;  (°)  —  neon  at  3  x  I0~*  T;  (A) 
»  krypton  at  5  x  IO'5  T;  and  (0)  -  xenon  at  3  x  10'5  T. 


unity  are  required  to  obtain  the  most  efficient  electron- 
beam  transport.  This  limits  the  maximum  instantaneous 
current  transport  efficiency  to  76%  and  the  charge  trans¬ 
port  efficiency  obtained  in  these  experiments  to  approxi¬ 
mately  56%  within  the  5.08-cm  diam  of  the  Faraday  cup. 
The  charge  transport  efficiency  does  not  appear  to  be  a 
strong  function  of  the  channel  gas  mass  for  the  present 
experimental  conditions. 

Transverse  beam  centroid  oscillations  have  been  ob¬ 
served  with  real  frequencies  and  qualitative  scaling,  with 
ft  and  the  channel  ion  mass  in  agreement  with  the  ex¬ 
pected  values  from  the  ion  hose  instability. 

The  data  indicates  that  the  REB-plasma  electron  two- 
stream  instability  can  become  dominant  for  measured  f, 
values  as  low  as  1.1  for  our  experimental  conditions. 
Transported  current  signals  (Faraday  cup)  show  that  the 
beam  current  also  begins  to  become  strongly  modulated 
when  ft»  1.  This  is  accompanied  by  a  loss  in  the  trans¬ 
port  efficiency  and  the  transported  beam  energy.  Both  the 
charge  transport  efficiency  and  transported  beam  energy 
peak  for  f,  values  around  unity  and  begin  to  fall  off  as/, 
is  increased  further.  The  initiation  of  the  two-stream  in¬ 
stability  for /,  values  slightly  less  than  unity  is  likely  due 
to  beam  electrons  interacting  with  plasma  electrons  on  the 
periphery  of  the  beam  radius. 
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The  results  of  this  study  can  be  compared  to  our  pre¬ 
vious  dkperiments  concerning  electron-beam  transport  in 
laser-induced  IFR  channels  utilizing  the  MELBA  elec¬ 
tron-beam  generator  [4]  at  the  parameters:  voltage  = 
—0.8  MV;  injected  beam  currents  from  200  to  300  A;  and 
pulse  lengths  variable  from  0.4  to  1  /is.  In  those  experi¬ 
ments  [8]  the  background  gas  pressure  was  much  higher, 
so  electron  ynpact  ionization  .increased  the  channel  ion 
density  during  the  pulse.  The  present  experiments  have 
operated  at  background  gas  pressures  low  enough  that 
REB-induced  electron  impact  ionization  is  negligible  dur¬ 
ing  the  pulse.  Peak  electron-beam  current  transport  effi¬ 
ciencies  in  the  previous  laser  ion  channel  guiding  exper¬ 
iments  were  as  high  as  80%  for  comparable  pulse  lengths, 
compared  to  76%  in  the  present  experiments.  These  ex¬ 
periments  demonstrate  that  ion  channel  guidance  is  an  ef¬ 
fective  means  of  transporting  long-pulse,  high-current 
electron  beams.  However,  the  neutralization  fraction  fe 
must  be  kept  to  less  than  unity  in  order  to  avoid  streaming 
instabilities. 
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Argon  ion  excitation  by  relativistic  electrons: 

I.  Collision  cross  sections  and  deposition  efficiencies 

By  D.  B.  McGARRAH*  and  M.  L.  BRAKE 
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Calculations  of  the  electron  impact  excitation  cross  sections  and  deposition  efficiencies 
for  singly  ionized  argon  with  electrons  of  energies  up  to  and  including  relativistic  values  have 
been  made  using  the  first  Bom  approximation  and  the  generalized  oscillator  formalism 
Deposition  efficiencies  for  fast  electrons  were  generated  from  the  Peterson  and  Green 
integral  equation.  Cross  sections  and  efficiencies  were  produced  for  29  transitions  from 
the  ground  state  configuration  of  Aril  to  excited  energy  levels  with  (Ne)3s23p44s  and 
(Ne)3s23p43d  configurations  and  for  40  transitions  between  excited  energy  levels  from  4s 
and  3d  to  4p  orbitals.  Efficiencies  are  constant  for  electron  energies  above  1  keV  to  10  MeV. 
Electrons  ejected  from  inner  shells  contribute  up  to  12ft  of  the  efficiency  of  the  transition 
for  electrons  above  10  keV. 


1.  Introduction 

Although  relativistic  electron  beams  (REB)  are  frequently  used  to  pump  gas  lasers  (e.g.t 
Werner  et  al.  1977),  the  basic  understanding  of  how  the  fast  electrons  populate  excited 
states,  particularly  for  ions,  is  still  only  approximately  known.  We  have  recently  observed 
unusual  light  emission  from  REB  produced  plasmas  in  rare  gases  (Brake  et  al.  1986;  Brake 
&  Repetti  1988).  In  order  to  understand  these  studies  of  argon  ion  emission  produced  by 
a  relativistic  electron  beam  it  is  necessary  to  generate  a  complete  set  of  cross  sections  and 
deposition  efficiencies  for  singly  ionized  argon  experiencing  collisions  with  relativistic 
electrons. 

Peterson  and  Allen  (1971)  have  calculated  efficiencies  for  certain  transitions  in  argon  for 
neutral  atoms  experiencing  collisions  with  electrons  of  energy  no  more  than  10  keV.  They 
used  a  semiempirical  approach  based  on  the  Bom  approximation  following  the  develop¬ 
ment  of  Mott  and  Massey  (1965).  This  assumes  that  spatial  variations  are  not  of  interest 
and  uses  a  continuously  slowing  down  approximation  to  calculate  the  energy  degradation 
of  the  primary  electron.  More  recently,  Bretagne  et  al.  (1986)  investigated  collisional  ex¬ 
citation  of  neutral  argon  due  to  relativistic  electrons.  In  this  paper,  the  methods  of  Peter¬ 
son  and  Allen  (1971)  and  Bretagne  et  al.  (1986)  are  used  to  generate  the  relevant  cross 
sections  and  efficiencies  for  singly  ionized  argon  experiencing  collisions  with  electrons  rang¬ 
ing  in  energy  from  a  few  eV  to  10  MeV. 

2.  Calculation  of  electron  impact  excitation  cross  sections 

A  quantum  mechanical  theory  of  inelastic  collisions  of  fast  charged  particles  with  atoms 
or  molecules  was  developed  by  Bethe  (1932)  using  the  Bom  approximation.  The  rdativis- 
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tic  version  of  Bethe’s  differential  cross  section  in  generalized  oscillator  strength  formalism 
for  inelastic  collisions  of  electrons  with  atoms  taken  from  the  review  by  Inokuti  (1971)  }s 

doji  =  WW'(k‘/k) 

dQ  h*c*  {tf2  -  (W-  W)2/h2c2]2  (1) 

where  au  is  the  cross  section  for  excitation  from  level  /  to  level  k  and  0  is  the  solid  angle 
subtended  by  the  trajectory  of  the  particle  before  the  collision  with  that  of  the  trajectory 
of  the  incident  particle  after  the  collision.  W  (If")  is  defined  as  the  total  energy  of  the  in¬ 
cident  particle  before  (after)  the  collision,  and  k  ( k ')  as  the  wave  number  of  the  particle 
before  (after)  the  collision.  The  physical  constants  e,  h,  and  c  are  the  charge  of  an  electron, 
Planck's  constant,  and  the  speed  of  light  in  a  vacuum,  respectively.  The  relativistic  form 
factor  i?  is  given  by  Bretagne  et  al.  (1986): 


ll*<*)|J 


QFU{K) 

Eik 

QFn(K)  (1  —P2)fikEiic 


(tot)  *  1 


(*h02)«l 


where  for  a  relativistic  incident  particle,  Q  is  defined  as  follows: 

h2K2  —  Efc/c2 
V  2m 

and 


Urn  £•*<*)=/*  (4) 

and 

K  =  k-k'  (5) 

where  HK  is  the  amount  of  momentum  transferred  by  the  collision,  F*(AT)  is  the  general¬ 
ized  oscillator  strength  for  the  transition  from  level  i  to  level  it  and  is  a  function  of  momen¬ 
tum  transfer  of  the  collisional  partners,  /*  is  the  optical  oscillator  strength  and  is  constant 
for  a  given  transition,  Ey,  is  the  energy  difference  between  levels  r  and  k,  and  0  is  the 
relativistic  beta,  i.e.,  the  ratio  of  the  speed  of  the  particle  to  the  speed  of  light  in  a  vacuum, 
and  m  is  the  mass  of  the  incident  particle.  The  total  cross  section  for  collisional  excitation 
can  be  found  by  integrating  equation  (1). 

First  the  differential  cross  section  with  respect  to  solid  angle  Q  of  the  change  in  trajec¬ 
tory  of  the  incident  electron  is  transformed  to  a  differential  cross  section  with  respect  to 
the  scattering  angle  0.  Applying  the  cosine  law  to  the  vectors  in  equation  (S)  allows  trans¬ 
formation  of  f(0)dd  tof(K)d(K2).  Then  the  total  cross  section  for  collisional  excitation 
can  be  found  by  integrating  the  differential  cross  section  over  the  limits  of  possible  momen¬ 
tum  transfer  using  the  values  of  the  relativistic  form  factor  at  the  limits  of  low  and  high 
values  of  (Ka$)  as  indicated  in  equation  (2).  The  resulting  expression  is 


SrajR2  fit 
mc20 2  Em 


(6) 


where  R  is  the  Rydberg  energy,  Oo  is  the  Bohr  radius,  and  C*  is  a  function  of  the  gener¬ 
alized  oscillator  strengths  which  in  practice  is  unknown  but  can  be  evaluated  by  experi¬ 
ments.  C*  is  generally  expressed  as  of  where  f  =  R/Em.  Bretagne  et  al.  (1986)  found  that 
a  multiplier  of  o  *  1/4  gave  good  agreement  with  experimentally  measured  values  for  neu¬ 
tral  argon.  It  is  assumed  that  this  will  hold  for  argon  ions  as  well  (see  “Note  added  in 
proof"). 
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Since  the  Born  approximation  is  only  valid  at  high  energies,  the  expression  for  the  ex¬ 
citation  cross  section  must  be  modified  so  that  it  can  be  used  at  lower  energies.  The  cross 
section  is  modified  by  defining  a,k  equal  to  Bik  times  the  value  of  a,k  given  in  equation  (6) 
where  Bik  is  given  as  (Bretagne  et  al.  1986): 


B,k  =  1  - 


2  Eik 
mc2&2  * 


(7) 


To  evaluate  the  excitation  cross  section  it  was  necessary  to  obtain  values  for  the  energies 
Em  and  the  oscillator  dipole  strengths/*.  Although  the  energy  levels  of  Aril  have  been 
well  documented  (Moore  1949),  it  is  more  difficult  to  find  extensive  tables  of  line  strengths 
for  many  transitions  in  ions.  Although  a  few  line  strengths  are  known  (Wiese  et  al.  1969), 
a  self-consistent  set  of  oscillator  line  strengths  for  all  transitions  of  interest  had  to  be  gen¬ 
erated  from  an  atomic  spectra  simulation  developed  by  Cowan  (1970).  The  Cowan  code 
consists  of  a  series  of  three  programs.  The  first  generates  single  configuration  radial  wave- 
functions,  the  second  is  an  interface  program,  and  the  last  program  performs  energy  plane 
wave  Born  collision-strength  calculations.  This  well-known  atomic  physics  model  can  be 
used  to  produce  energy  levels  and  line  strengths  for  all  allowed  transitions  between  all  en¬ 
ergy  levels  of  all  given  electron  configurations  from  one  parity  to  the  opposite  parity. 

For  convenience  the  cross  section,  including  the  modifier  from  equation  (7),  will  be  re¬ 
written  in  the  following  form: 


(8) 

,  2.55  X  10-19  ,  j 

A,  =  [cm’eVl, 

(9) 

9.96  X  10"“  r  2„ 

A2  = - - [cm2!2, 

(10) 

A 3  =  1.02  x  10«teVla 

(ID 

Ak  =  0J. 

(12) 

The  parameters  Eik  and/*  as  determined  from  the  Cowan  code,  are  given  in  table  1  and 
table  2  for  all  the  transitions  that  were  examined.  Measured  oscillator  strengths  compiled 
by  Wiese  et  al.  (1969)  are  also  given  in  table  1  and  table  2.  Taking  into  account  the  fact  that 
the  measured  oscillator  strengths  have  an  error  of  about  30-50%,  there  is  excellent  agree¬ 
ment  between  the  measured  values  and  those  generated  by  the  Cowan  code. 

The  cross  sections  presented  in  figures  1-4  are  calculated  using  equation  (6)  multiplied 
by  fi*  as  given  in  equation  (7)  to  allow  for  consideration  of  low  energy  as  well  as  relativis¬ 
tic  electrons.  Note  that  the  cross  sections  in  figures  1  and  2  are  for  transitions  originating 
from  the  ground  configuration  of  the  argon  ion  and  those  in  figures  3  and  4  are  for  tran¬ 
sitions  beginning  from  the  (3s*3ps)  4s  configuration,  i.e.,  it  is  presumed  that  the  level  al¬ 
ready  exists.  Only  allowed  transitions  were  considered.  We  have  generated  cross  sections 
for  all  transitions  listed  in  table  1  and  table  2.  The  selection  rules  applied  were: 

AL  =  ±1  A/  =  0,  ±1  (13) 


where  L  is  the  magnitude  of  the  orbital  angular  momentum  of  the  ion  and  J  is  the  magni¬ 
tude  of  the  total  (orbital  plus  spin)  angular  momentum.  The  transitions  3p-4s  and  As-Ap 
were  chosen  to  be  displayed  as  they  are  considered  to  be  significant  transitions,  for  not  only 
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Table  I.  Energy  (£**)  end  oscillator  line  strengths  (/#)  for  transitions  in  singly  ionized 
argon  between  the  ground  levels  with  (3s23p5)  configuration  and  the  excited  levels 
with  configuration  (3s23p4)  4s,  id.  Whenever  possible  the  line  strengths 
determined  from  the  Cowan  code  were  compared  with  experimental 
measurements  (Wiese  et  at.  1969). 


Transition 

Elk 

U 

Cowan 

fik 

(NBS) 

1. 

3pl2P3n\ 

-As  [*P3n\ 

16.43 

2.2Se-4 

2. 

3p[2P3/2) 

~As  l*Pi/2] 

16.53 

3.35«-3 

2.5e-3 

3. 

3p[2P3/2] 

-As  I 'Pm] 

16.59 

2-5<-5 

1.5e-4 

4. 

3p[2P3/2] 

-As  [2P3/2] 

16.91 

0.1920 

0.18 

5. 

3p[2Pi/2\ 

-As  l2Pin] 

17.04 

0.0404 

0.037 

6. 

3p[2P3/2\ 

-4s'(2£>3/2] 

18.30 

3.78e-2 

7. 

3p[2P3/2'\ 

-As'[2Di/2] 

18.33 

0.0004 

8. 

3pl^3/2l 

-4s'[2S1/2] 

21.16 

0.0697 

9. 

3p[2P3n] 

-3d  l2D3/2] 

18.72 

0.0738 

10. 

3pi2P3n\ 

-3dl2Di/2) 

18.78 

0.1527 

11. 

3p[2Pin] 

-  3d'l2Ds/2] 

21.50 

0.5483 

12. 

3p[2P3n\ 

-3d,[2D3/2] 

21.55 

0.IS83 

13. 

3pl2P3/2] 

-3d'[2Pi/2] 

21.65 

0.2870 

14. 

3p[2P3/2) 

-3d‘[2Pin] 

21.68 

0.4830 

15. 

3pl2P3/2) 

-3^1%^] 

22.83 

0.4830 

16. 

3pl2P3,2] 

-3d'[2Di/2) 

22.54 

0.0550 

17. 

3pl2P3/2) 

-Zd'frhni 

22.56 

1.15e-3 

18. 

3pl2Pu2] 

-As  [*P3/2] 

16.35 

0.0007 

9.9e-4 

19. 

3pl2Pt/2 1 

-As  l*Pl/2] 

16.42 

8.5e-4 

6.1e-4 

20. 

3pl2Pl/2) 

-As  l2P3/2] 

16.74 

0.0642 

0.072 

21. 

3p[2Pin\ 

-As  llPm  1 

16.86 

0.1464 

0.15 

22. 

3pl2Pxn\ 

-As'lllhn] 

20.98 

0.0740 

24. 

3p[2Pini 

-4j'[2Si/2] 

18.54 

0.1883 

25. 

3pl2P3/2\ 

-3d'[2Z)va] 

21.38 

0.3856 

26. 

3pl2Pt/i] 

-  3d'l2Pi/2) 

21.48 

0.3238 

27. 

3pl2Pl/2) 

-3dl2Pt/2] 

21.51 

0.2893 

28. 

3pl2P3/2) 

-3d'l2Sl/2] 

22.81 

0.5290 

29. 

3pl2Pt/2] 

—  3d,l2D3/2] 

22.38 

9.465 e-2 

the  experimental  studies  mentioned  earlier  (Brake  et  al.  1986;  Brake  &  Repetti  1988),  but 
for  the  general  case  of  transitions  in  argon  ion  lasers. 

It  is  difficult  to  compare  the  calculated  cross  sections  to  experimental  values  since  very 
little  data  exists  for  excited  states  of  argon  ions,  particularly  at  relativistic  energies.  How¬ 
ever,  a  comparison  between  the  calculated  values  of  the  excitation  cross  section  and  those 
measured  by  Feltsan  and  Povch  (1970)  for  the  4 s’-4p'  transition  (4277  A)  upon  radiative 
decay  from  the  upper  to  lower  level  at  low  energies  is  made  and  shown  in  figure  5.  Unlike 
experiments  examining  excitation  of  neutral  argon,  it  is  not  possible  to  directly  compare 
excitation  cross  sections  to  experimental  measurements  for  ions  since  the  calculation  pre¬ 
sumes  an  initial  target  of  Aril  already  in  a  particular  energy  level  which  is  generally  not  the 
case  experimentally  (most  experiments  begin  with  ground  state  neutral  argon  atoms).  Felt¬ 
san  and  Povch  (1970)  conducted  their  experiment  by  injecting  electrons  into  neutral  argon 
in  order  to  observe  transitions  in  singly  ionized  argon,  that  is,  they  assume  a  direct,  one 
step  ionization  and  excitation.  Our  calculations,  on  the  other  hand,  involve  a  three  step  ex¬ 
citation  process  as  follows: 
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Table  2.  Energy  (£,*)  and  oscillator  line  strengths  (/*)  for  transitions  in  singly  ionized 
argon  between  the  excited  levels  with  (3s23p5)  4s  and  3d  configuration  and  the 
excited  levels  with  configuration  (3 s*3p4)  4p.  Whenever  possible  the  tine 
strengths  determined  from  the  Cowan  code  were  compared  with 
experimental  measurements  (Wise  et  at.  1969). 


Transition 

Eik 

fit 

Cowan 

fit 

(NBS) 

1. 

4s  [*Pin\  -4 p  [*Pi/2] 

2.521 

0.1502 

.13 

2. 

4s  l*P3,2]  -4p  [4Din] 

2.846 

0.0228 

0.0165 

3. 

4s  l4P3n)  —  4p  l*Pi/2) 

2.419 

0.0546 

0.0254 

4. 

4s  [4Pm]  —  4p  [*Din\ 

2.782 

0.3439 

0.299 

5. 

4s  1<P1/21  —  4p  (<P3/21 

2.355 

0.2024 

0.170 

6. 

4s  1 2P3/2]  -4 p  f2P1/2] 

2.633 

0.1517 

0.132 

7. 

4s  l2P3/2]  —  4p  [2P3/2] 

2.702 

0.1610 

0.128 

8. 

4s  l2P,/zl  -4p'[2Pin] 

4.350 

0.0679 

9. 

4s  [2P1/2]  —  4p  I2P,/2] 

2.511 

0.0661 

0.057 

10. 

4s  [  2P1/2)  —  4p  [2P3/2] 

2.579 

0.4004 

0.391 

11. 

4s  [2P1/2]  -4p'(2P3/2) 

4221 

0.0270 

12. 

4s'[  2D3/2)  -4p  I  2Pj/2] 

1.247 

0.0279 

13. 

4s'(2D3/2l-4p'l2£»3/2l 

3.099 

0.0203 

14. 

4s'  ( 2D2/2l  -4p'[2P3/2] 

2.964 

0.0280 

0.057 

15. 

4s'I2Dj/2]  —  4p  [2P3/2] 

1.281 

0.0363 

16. 

4s'[2Dj/2]  -4p'[2Z>J/2] 

3.064 

0.1455 

0.142 

17. 

4s'{2Z>j/z]  -4p'{2P7/2] 

2.735 

0.3968 

0.387 

18. 

4s'l2DJ/2]  -4p'[2P3/2l 

2.928 

0.1499 

0.18 

19. 

4s'12S1/21  —  4p'l2Pj/2] 

1.805 

0.2339 

0.37 

20. 

4s'[2Si/2J  -4p'I2P3/2) 

1.812 

0.4693 

0.75 

21. 

3d[2D3/2\-4p'[2Ds/2\ 

2.680 

0.0357 

22. 

3d[2D3/2]-4p'[2Pi/2\ 

2.545 

0.0336 

23. 

3dl2Ds/2]  -4p'f2Z)J/2) 

2.619 

0.2615 

24. 

3rfI2D3/2]  -4p'[2F7/2) 

2.291 

0.1210 

25. 

3<f  C2Os/2l-4p'C2P3/2] 

2.484 

0.0372 

26. 

3d'1 2£>3/2J  -4p  I2P3/2] 

1.884 

0.1508 

27. 

3<f'l2Z)S/2l~4p'(2P3/2l 

1.903 

0.0175 

28. 

3<f'(^/zl  -4p  l\Pl/2 1 

2.008 

0.1044 

29. 

3rf'l2£V2i  ~4p  l2P3/2l 

1.939 

0.0104 

30. 

3<#'(2Pvj]  -4p  [2P3/2] 

2.042 

0.0471 

31. 

3rf'f2Pj/zl  -4p'(2P3/2J 

0.3942 

0.0244 

32. 

3d'(2P,/zl  -4p  (2Pl/2l 

2.138 

0.1638 

33. 

3<f'(2P,/zl  -4p'(2P3/2l 

0.4216 

0.0111 

34. 

3<n2S,/zi  "4p  (2P3/2J 

8.450 

0.0148 

35. 

3d'[2S1/2)  -4p'[2P3/2] 

6.820 

0.6865 

36. 

3d'[2S1/2]  -4 p'[2P3/2) 

5.094 

0.0108 

37. 

3d'l2Ds/2)  -4 p  l2P3/2) 

2.925 

0.0626 

38. 

3dVpsn\  -4p'l2P3/2l 

0.4305 

0.0327 

39. 

3d'(2D3/2l-4p  (2Pl/aJ 

3.007 

0.0447 

40. 

3d't2DJ/2]  -4p'I2P,/jl 

0.4122 

0.0266 

e  +  Arl  -  ArII(3p)  +  e  (14) 

e  +  ArII(3p)  -»  Aril  (4s)  +  e  (15) 

e  +  Aril  (45)  -»  ArII(4p)  +  e  (16) 

followed  by  radiative  decay 

ArII(4/7)  -  ArII(4j)  +  hr.  (17) 

?-< 
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Figure  1.  Excitation  cross  sections  as  a  function  of  electron  energy  for  transitions  from  the  / = 3/2 
level  of  the  ground  configuration  (3s23p5)  of  Aril  to  the  (3s23p4)  4s  configuration  where 
the  specific  transitions  are  as  foUows:  3pl2Pi/2]  -  4sl2P3/2)  (A);  3p[2P3/2]  -  4s[2PI/2J  (B); 
3pl2P3/2)  -4s[4P3/2]  (Q;  3pf 2P3/i\  —  4j(4P$/jJ  (D);  and  3p{2PJ/2)  -4 s[4PI/2]  (E). 


These  steps  involve  the  calculation  of  the  probability  of  ionization  and  the  probability  of 
a  transition  from  level  i  to  level  k  given  a  specific  /  and  k  being  one  of  the  allowed  final 
levels  as  well  as  the  branching  ratios  for  reaction  (17).  These  various  probabilities  have  been 
incorporated  into  the  calculated  effective  cross  sections  shown  in  figure  5. 


Figure  2.  Excitation  cross  sections  as  a  function  of  electron  energy  for  transitions  from  the  7*  1/2 
level  of  the  ground  configuration  (3s23p!)  of  Aril  to  the  (3s23p4)4s  configuration,  where 
the  specific  transitions  are  as  follows:  3p[2P\n\  -  4s{2Pi/j1  (F);  3p[ 2P\n\  -  4s[2P3/2J  (G); 
3pl2/»,/jl  -4sI4Pl/2J  (H);and3pl2Pl/2l  -4j(4P,/IJ  (I). 
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Figure  3.  Excitation  cross  sections  as  a  function  of  electron  energy  for  transitions  in  Aril  from 
(3s23ps)4s(  4PJ  levels  to  (Sr^p^pl4!),4/*]  levels  where  the  specific  transitions  are  as  fol¬ 
lows:  4s(  4PI/2J  -4pI401/2J  (J);4j{4Pj/j1  -4pf4Z>,/jl  (K); 4f[ 4/>s/2]  ~4p[*Pin)  (L);  4r[  4PJ/2] 
-4p[*Pin\  (M);  and  4r{4PI/2]  -4pl4/>3/2]  (N). 

As  an  additional  check,  figure  6  shows  the  low  energy  calculated  cross  section  for  the 
4s'-4p’  transition,  i.e.,  equation  (16),  compared  to  the  common  general  expression  for  elec¬ 
tron  excitation  given  in  plasma  formularies  (Book  1983): 


Oik 


ik  =  2.36  x  10~u  — — —  lcm2l 


(18) 


Figure  4.  Excitation  cross  sections  as  a  function  of  electron  energy  for  transitions  in  Aril  from 
<3r23p5)4r[2P]  levels  to  (3s23p5)4pI2P)  levels,  where  the  specific  transitions  are  as  follows: 
4r[  2P1/a)  - 4p(  2Pin)  (O);  4s(  2Pin)  - 4p{  2Pin]  (P);  4r{  lPln]  -  4pl  *Pl/2)  (Q);  and 4r{ 

~  4p(  2P|/j]  (R). 
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Figure  5.  Effective  excitation  cross  section  for  three  step  process  of  equations  (14H16)  as  a  func¬ 
tion  of  electron  energy  for  the  4s'  -  4p'  transition  in  singly  ionized  argon  (4277  A).  The  (•)  points 
indicate  the  cross  section  values  measured  experimentally  by  Fehsan  and  Povch  (1970). 
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The  Petersor. 
Down  Approach 
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\ 


where and  Em  are  again  the  line  strength  and  energy  of  the  transition  and  here  <  is  the 
energy  of  the  electron  and  g(i,k )  is  the  Gaunt  factor  (generally  —0.2  for  ions)  (Ono  1980). 
(Cross  sections  for  other  transitions  were  also  compared  to  their  expression  and  indicated 
similar  agreement.)  Note  that  this  approximate  expression  is  inversely  proportional  to  en¬ 
ergy  and  continues  to  increase  to  infinity  as  energy  goes  to  zero.  Physically  a  transition  can¬ 
not  occur  for  values  of  energy  transfer  less  than  the  difference  in  energy  between  the  upper 
and  lower  levels  of  the  transition.  Equation  (18)  thus  leads  to  large  errors  at  values  of 


Fiouu  6.  Excitation  cross  section  as  a  fumsion  of  electron  energy  for  the  4s'-4pb»aw^)iti0(1 

culated  from  equation  (8)  (smooth  line)  compared  with  the  cross  sections  generated  from 
(18)  (dashed  line).  Note  that  the  common  expression  (dashed  line)  fails  at  energies  near  u* 
old  for  the  transition. 
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energy  transfer  which  are  near  the  threshold  for  the  transition  as  well  as  energy  levels  in 
the  relativistic  regime. 


3.  Calculation  of  electron  beam  deposition  efficiencies 


The  Peterson  and  Green  (1968)  integral  equation  based  upon  the  Continuously  Slowing 
Down  Approximation  (CSDA)  was  used  to  calculate  the  energy  apportionment.  In  this  way 
the  number  of  excitations  generated  by  the  beam  electron  and  all  subsequent  generations 
of  secondary  electrons  can  be  counted.  The  number  of  transitions  from  energy  level  /  to 
k  produced  by  a  single  beam  electron  with  kinetic  energy  T  is  given  by  (Peterson  &  Green 
1968): 

/•(T-/y)/ 2 

N*(T)  =N,UT)  +  2  Na{X)nu(T,X)dX  (19) 

J  •'£,* 


with 


Q,*(X)dX 

LAX) 


(20) 


and 


nu(T,Ts) 


fT  SU(X,  T,)dX 
'27i+/y  LAX) 


(21) 


where  Nu(T)  is  equal  to  the  number  of  transitions  i  to  k  produced  by  the  beam  electron 
and  all  subsequent  generations  of  electrons  generated  and  N*(T)  is  equal  to  the  number 
of  transitions  /  to  k  produced  directly  by  a  primary  electron  with  kinetic  energy  T.  The 
lower  limit  is  the  minimum  amount  of  energy  that  must  be  transferred  to  the  ion  to  cause 
a  transition  from  level  /  to  level  k.  The  upper  limit  is  the  kinetic  energy  of  the  incident  elec¬ 
tron.  The  term  nu(T,  Ts)  is  equal  to  the  number  of  secondaries  per  unit  energy  with  kinetic 
energy  Tt  generated  from  ionization  of  theyth  shell  of  the  atom/ion  in  energy  level  i  by  an 
electron  with  kinetic  energy  T.  The  lower  limit  is  the  minimum  amount  of  energy  which 
must  be  transferred  from  the  incident  electron  with  kinetic  energy  T  to  the  ion  so  that  an 
electron  with  kinetic  energy  T,  will  be  ejected.  Again,  the  upper  limit  is  the  kinetic  energy 
of  the  incident  electron. 

The  total  loss  function  for  an  electron  experiencing  inelastic  scattering  off  of  the  target 
atom/ion  in  energy  i  level  can  be  written  as  (Bretagne  el  at.  1986): 


/•(7>ztf)/2  da  IT  El 

LAT)  =  ZE«**(r)  +  ZBj  E  dE.  (22) 

k  J  Jl(j 

It  includes  both  excitation  and  ionization.  Again,  modifications  have  been  made  to  account 
for  low  energy  electrons  with  By  defined  as  follows: 


By  =  1 


2/y 

~  mc202 


(23) 


where  Iu  is  the  ionization  potential  for  an  atom/ion  in  energy  level  i  to  lose  an  electron  in 
shell  j.  Iu  is  defined  in  this  manner  so  that,  in  the  valence  shell  for  an  excited  level,  lu  is 
equal  to  the  ionization  potential  minus  the  energy  difference  between  the  ground  state  and 
that  level.  For  inner  shells,  the  ionization  potential  remains  the  same  since  all  excitation 
transitions  under  consideration  involve  valence  electrons  only.  Sy(T,E)  is  the  differential 
(in  energy)  ionization  cross  section  for  an  incident  electron  with  kinetic  energy  T  transfer- 
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ring  energy  E  to  the  target  atom/ion  in  energy  level  i  by  ejecting  an  electron  from  the 
yth  shell.  It  was  calculated  as  the  sum  of  two  terms  representing  glancing  and  hard  colli¬ 
sions.  The  glancing  collision  term  is  an  extension  of  the  differential  excitation  cross  section 
where  the  excitation  energies  are  allowed  for  transitions  into  continua  and  is  recast  from 
f(K)d(K 2)  to  the  form/(Q)d(In  Q).  The  glancing  component  of  the  differential  ioniza¬ 
tion  cross  section  is  thus: 


•0/  \nj(K,T)\2d\]nQ\ 

Q 


(24) 


„  me2/ 3 

Evaluating  the  integral  and  using  the  result  of  Bretagne  et  al.  (1986)  for  the  lower  limit 

T2(l  -02) 


e.= 


2mcz/33 


(25) 


The  upper  limit  on  Q,QJt  for  a  glancing  ionization  collision  ejecting  an  electron  in  the  yth 
shell,  is  given  by  Bretagne  et  al.  (1986)  for  neutral  argon.  They  show  that  Qj  —  1500  eV  for 
the/T shell,  lOOeV  for  the  L  shell,  and  1.1  +0.3(1.  -exp(-(£'-/tt#)/5)]  for  the  Af  shell 
where  A#  =  A#  —  £<>  A#  equal  to  27.63  eV  is  the  ionization  potential  of  Aril  in  the  Af 


Figure  7.  Deposition  efficiencies  of  the  <  3sz3p4)  4 s  and  3 d  levels  originating  from  the  J  =  3/2  level 
of  the  ground  configuration  (3sz3p5)  of  Aril  as  a  fraction  of  the  incident  energy  where  the  final 
states  are  as  follows:  4sIz/>j,21  (A);  4j(4/>j/2]  (B);  3d'l2D5/2]  (C);  3d'(2SJ/2J  (D);  id'^Pj/i) 
(E);  3d'(2P|/jl  (F);  3d[2DJ/2]  (G);  4s(4P5/2]  (H);  4s(zP,/jl  (I);  3dT2Oj/jl  (J):  ts'^Dwi) 
<K);  3d(  2D1/2)  (L);  4r{4Pl/2]  (M). 
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shell,  and  Et  is  the  fth  energy  level  of  the  ion.  The  unsubscripted  E  refers  to  the  energy 
transferred  from  the  incident  electron  during  the  collision.  We  assume  that  these  limits  hold 
for  argon  ions,  as  well  as  neutrals,  if  the  appropriate  ionization  potential  is  used  for  the 
M  shell  of  the  ion.  The  glancing  component  of  the  differential  ionization  cross  section  is 
then  given  by: 


dap  _ 
dE  t 


(26) 


The  hard  collision  is  considered  to  be  a  scattering  of  the  fast  electron  by  a  free  electron 
at  rest.  The  expression  for  this  differential  cross  section  has  been  calculated  by  Bretagne 
etal.  (1986): 
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Figure  8.  Deposition  efficiencies  of  the  (3s*3p4)  4s  and  3d  levels  originating  from  the  7=  1/2  level 
of  the  ground  configuration  <3s23ps)  of  Aril  as  a  fraction  of  the  incident  energy  where  the  final 
states  are  as  follows:  4j(2P3/2]  (N);  4*llPl/a]  (O);  3d'[2S3/2]  (P);  3d'[2Pl/2)  (Q);  3d[  iDin\ 
(R);  4j'12Dj/2]  (S);  3d'[ 2P3/j]  CD:  4sl4/*,^]  (U);  4il4P3/2l  (V);  3d'[JfVal  <W). 
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where  Nj  is  the  number  of  bound  electrons  in  they  th  shell  and 

£&y  _  °pj 
dE  4*aaoR 


(29) 


where  <rw  is  the  photoionization  cross  section  of  theyth  shell  and  a  is  the  fine  structure  co¬ 
efficient.  The  term  dfj/dE  goes  as  E~i  S  for  large  values  of  E  (see  Bretagne  et  at.  1986).  It 
is  assumed  that  the  inner  shells  of  singly  ionized  argon  have  the  same  energy  dependence 
as  those  of  neutral  argon.  Since  dfj/dE  decreases  with  increasing  energy,  an  upper  limit  on 
its  value  can  be  set  by  calculating  its  value  at  the  lower  limit  of  energy,  the  ionization  po¬ 
tential.  This  can  be  determined  by  using  equation  (29)  and  the  value  of  the  photoioniza- 
tion  cross  section  at  the  ionization  potential  as  given  by  Marr  (1968). 

The  beam  deposition  efficiencies  are  defined  as  the  fraction  of  the  energy  of  the  electron 
beam  deposited  in  the  gas  by  causing  transitions  from  energy  level  /  to  k: 


V,*(T)  m 


Nik(T)Eik 

T 


(30) 


Efficiencies  are  very  useful  in  chemical  kinetic  models  of  REB  produced  plasmas  and  lasers. 
In  particular,  a  pumping  rate  coefficient  for  raising  ions  in  energy  level  f  to  a  higher  en¬ 
ergy  level  k  due  to  collisions  with  relativistic  electrons  can  be  defined  as 
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Figure  9.  Deposition  efficiencies  of  the  <3s23p4)  4 p  levels  originating  from  the  /  *  3/2  level  of 
the  ground  configuration  (3r23p5)  of  Aril  as  a  fraction  of  the  incident  energy  where  the  final 
states  are  as  follows:  4 p{  2P|/j]  (A);  4pl  2Pj/j]  (B);  4p'(  2Pj/i)  (C);  *P l  *Di/i)  (W;  4pt  P3  2I 
(E);  4p'(2Pj/j]  (F);  4p'(  2Dj/j]  (O);  4p'|2Pl/2l  (H);  4p'(lP7/2l  (D 
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where  Jb  is  the  current  density  of  the  electron  beam,  L,{T)  is  the  loss  function  for  an  elec¬ 
tron  with  kinetic  energy,  T,  incident  on  a  target  in  energy  level  i,  and  £,*  is  the  energy  of 
the  transition. 

Efficiencies  for  transition  in  Aril  are  displayed  in  figures  7-10.  For  most  transitions  the 
values  of  the  efficiencies  are  roughly  constant  above  a  beam  energy  of  1  keV,  although  in 
a  few  cases  ij  continues  to  increase  up  to  beam  electron  energies  of  10  keV.  The  results  of 
these  efficiency  calculations  will  be  used  to  model  the  plasma  chemistry  of  argon  gas  ex¬ 
cited  and  ionized  by  relativistic  electrons  (McGarrah  &  Brake  1990). 

The  efficiency  for  ionization  can  be  considered  by  examining  the  extreme  case  of  level 
/  equal  to  the  ground  level  and  level  k  going  to  infinity.  It  would  have  a  form  analogous 
to  that  for  excitation  given  in  equation  (30).  Since  the  energy  per  ion  pair  is  given  by  {T/N), 
which  is  equal  to  (Eik/ri(T)),  where,  for  this  case,  Nis  the  number  of  transitions  from 
ground  to  infinity,  i.e.,  the  number  of  ionizations,  and  £*  is  simply  the  ionization  poten¬ 
tial.  Energy  per  ion  pair  values  for  neutral  argon  are  known  to  be  fairly  constant  (Bekefi 
1976).  Since  the  ionization  potential  is  constant,  it  is  reasonable  that  the  efficiency  i j(7”) 
should  also  be  roughly  constant.  Ejection  of  inner  shell  electrons  contributes  approximately 
from  5%  to  as  much  as  12<7o  of  the  transitions  for  electron  energies  higher  than  10  keV.  A 
few  exceptions  have  lower  contributions  for  inner  shell  electrons. 


Fioujue  id.  Deposition  efficiencies  of  the  (3s23p4)  4 p  levels  originating  from  the /=  1/2  level  of 
the  ground  configuration  (3s13pi)  of  Aril  as  a  fraction  of  the  incident  energy  where  the  final  states 
are  as  follows:  4plJPJ/2]  (J);  4p\2Pxn\  (K);  4p  [2Pin]  <L);  4p\*Din)  <M);  4p’{1Pin\  (N); 
4p'l2Dj/jl  (O);  4p’[2Pxn\  (P);  4p[*Fin\  (Q). 
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4.  Conclusions 

A  complete  set  of  excitation  cross  sections  has  been  generated  for  69  allowed  transitions 
in  singly  ionized  argon.  The  generalized  oscillator  strength  model  for  relativistic  electrons 
was  used.  A  modifier  allowed  for  extension  to  lower  energy  electrons.  A  self-consistent  set 
of  oscillator  dipole  line  strengths  was  generated  for  Aril  by  the  Cowan  code.  When  appro¬ 
priate  corrections  were  made  it  was  possible  to  compare  calculated  cross  sections  with  those 
measured  in  experiments  starting  with  neutral  argon.  This  analysis  was  performed  for  a 
As'-Ap'  transition  measured  by  Feltsan  and  Povch  (1970)  and  compared  well  with  the  ex¬ 
perimental  values.  Efficiencies  for  populating  excited  states  were  calculated  by  using  the 
method  of  Peterson  and  Green  (1968)  to  quantify  the  effect  of  a  beam  electron  and  all  sub¬ 
sequent  generations  of  electrons  produced.  Results  in  most  cases  indicate  that  the  efficiency 
is  roughly  constant  above  1  keV.  Inner  shell  ionization  accounts  for  about  5  to  12%  of  the 
efficiency  when  considering  electrons  with  energies  greater  than  10  keV. 
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Note  added  la  proof 

Varying  a  from  0. 1  to  1 .0  for  an  energy  of  2  eV,  typical  of  observed  transitions,  the  term 
ln(2 mc2Cik/R)  changes  from  10.8  to  8.5  so  the  calculation  is  not  sensitive  to  the  value  of  a. 
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A  model  of  an  electron  beam  interacting  with  neutral  argon  was  developed  by  solving 
the  chemical  kinetic  rate  equations  for  the  time-dependent  populations  of  ground  and 
excited  levels  of  Arl  and  Aril  as  well  as  the  populations  of  electron  energy  groups.  In¬ 
tensities  of  spectral  lines  were  calculated,  from  predicted  population  densities  and  Ein¬ 
stein  coefficients,  and  compared  to  experimental  results. 

The  thermal  plasma  is  generated  during  the  beam  pulse  and  persists  for  some  time  after 
the  pulse  is  terminated.  Low  energy  levels  of  Aril  with  a  4s  or  Ap  valence  electron  have 
similar  time  profiles  to  the  plasma  density.  However,  high  energy  levels  of  Aril  with  a 
4 P"  or  4 p'  valence  electron  have  similar  temporal  profiles  as  the  beam  current.  They  are 
populated  predominantly  by  direct  ionization  and  exceed  levels  populated  by  thermal  elec¬ 
trons.  The  Ap’  levels  decay  rapidly  to  4s'  levels,  which  have  very  high  beam  deposition  ef¬ 
ficiencies  for  being  pumped  back  to  the  4p*  levels.  In  this  way,  certain  transitions  such  as 
4s’ -Ap'  in  Aril  are  pumped  directly  by  the  electron  beam.  Without  direct  ionization,  these 
transitions  would  be  negligible  compared  to  those  transitions  which  have  lower  energy 
levels. 


1.  Introduction 

Plasmas  produced  by  relativistic  electron  beams  (REB)  often  emit  large  amounts  of  light 
(Wei  el  dl.  1977).  This  property  has  allowed  them  to  be  successfully  used  as  pump  sources 
for  gas  lasers,  particularly  those  utilizing  atomic  transitions  in  rare  gases.  Experiments  per¬ 
formed  at  the  Intense  Beam  Interaction  Laboratory  of  the  University  of  Michigan,  search¬ 
ing  for  possible  lasing  lines  in  argon  ions,  observed  unusual  spectral  emissions  (Brake 
et  al.  1986;  Brake  &  Repetti  1988)  that  corresponded  to  transitions  in  Aril.  The  emissions 
that  were  unusually  intense  corresponded  to  transitions  originating  from  high  energy  lev¬ 
els  of  Aril  and  indicated  a  larger  population  of  certain  energetic  excited  levels  of  Aril  than 
would  be  indicated  assuming  a  thermal  distribution  of  excited  level  populations.  A  kinetic 
model  of  the  plasma  chemistry  was  needed  that  would  include  the  effect  of  beam  and  cas¬ 
cade  electrons  as  well  as  the  thermal  plasma  and  would  not  depend  on  assumption  of  par¬ 
ticular  equilibrium  conditions  in  order  to  understand  these  unusually  bright  spectral  lines. 


2.  Model 

In  the  first  paper  of  this  series  (McGarrah  &  Brake  1990),  an  “exact"  calculation  of  the 
energy  deposition  from  the  beam  was  made  in  the  sense  that  no  thermal  equilibrium  con¬ 
ditions  were  required.  As  discussed  in  that  paper,  a  Continuously  Slowing  Down  Approx- 

•D.B.  McGarrah’!  current  address  to  Sandia  National  Laboratories,  Division  8233,  Livermore,  CA. 
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imation  (CSDA)  was  utilized  for  calculation  of  the  energy  deposited  in  the  gas  due  to  the 
degradation  of  the  beam  and  subsequent  cascade  electrons. 

In  this  paper,  the  rate  equations  for  the  populations  of  Arl,  Arl*,  Arl(meta),  Aril, 
Aril*,  ebt  and  will  be  derived,  where 


Arl 

=  ground  level  neutral  argon. 

Aril 

=  ground  level  singly  ionized  argon, 

Arl* 

=  excited  neutral  argon, 

Arl(meta) 

=  metastable  neutral  argon, 

Aril* 

=  excited  singly  ionized  argon, 

eb 

=  beam  electrons, 

=  secondary  cascade  electrons. 

Although  optically  forbidden  transitions  generally  occur  with  much  lower  probability 
than  optically  allowed  transitions,  there  are  certain  metastable  levels  of  Arl  that  are  rela¬ 
tively  long  lived  and  can  subsequently  be  ionized  to  form  Aril.  The  population  of  metasta¬ 
ble  argon  is  expected  to  be  quite  small  in  comparison  to  that  of  the  ground  level  neutral 
argon,  but  for  completeness  it  will  be  included  as  a  possible  mechanism  for  producing  the 
singly  ionized  argon  (McGarrah  1989). 

The  spectral  emissions  of  Aril  observed  in  experiments  (Brake  et  ol.  1986;  Brake  & 
Repetti  1988)  are  identified  as  those  corresponding  to  the  transitions  (Ne)  3s*3p44s-»  (Ne) 
3s23p44p  and  (Ne)  3s23p43<f-*  (Ne)  3*23p44p.  Electron  configurations  with  a  principal 
quantum  number  above  4  were  not  inducted  since  no  transitions  from  higher  electron  con¬ 
figurations  were  experimentally  observed.  Only  optically  allowed  transitions  were  consid¬ 
ered.  The  selection  rules  applied  were: 

AL  =  ±1  AJ  —  0,±l  (1) 

where  L  is  the  magnitude  of  the  orbital  angular  momentum  of  the  ion  and  /is  the  magni¬ 
tude  of  the  total  (orbital  pulse  spin)  angular  momentum.  The  selection  rules  also  require 

A/=±l  (2) 

for  the  electron  in  transition  where  /  is  the  orbital  angular  momentum  for  the  single 
electron. 

The  ground  configuration  of  Aril  (ls22s23s23p5)  has  odd  parity,  i.e.,  the  sum  of  the 
angular  momentum  quantum  numbers  (/)  for  each  electron  is  equal  to  an  odd  number. 
Therefore  all  energy  levels  reached  from  the  first  excitation  must  have  even  parity  since  par¬ 
ity  is  required  to  change  during  a  transition.  The  energy  levels— in  Racah  notation  as  given 
in  Lengyel  (1966)— included  in  the  modd  are  given  in  table  1. 

Transitions  from  the  ground  configuration  to  (Ne)  3s23p44s  or  (Ne)  3s23p43<f  config¬ 
urations  with  oscillator  line  strengths  /#  less  than  0.01  are  not  used  in  the  modd  since  the 
stronger  lines  have  values  of  /*  on  the  order  of  1.0.  An  exception  to  this  is  made  only  if 
an  energy  level  is  the  lower  level  of  a  transition  observed  in  the  emission  spectra.  It  may 
have  a  high  transition  probability  of  being  reached  by  a  three  level  system  with  an  upper¬ 
most  levd  which  is  more  likdy  to  be  populated  through  a.transition  or  series  of  transitions 
from  the  ground  configuration  of  Aril  or  by  direct  excitation  during  the  ionization  of  Arl. 

The  physical  mechanisms  included  in  the  modd  for  the  populations  of  Arl  and  Aril  were 
the  following: 
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Table  1.  Energy  levels  of  Aril  included  in  the  model 


Configuration 


Term 


J 


a)  Odd  parity  energy  levels  with  ground  level  electron  configuration 
3r23p5  •  1 2P]  3/2.  1/2 


b)  Even  parity  energy  levels  with  excited  electron  configurations 


3*'3  p6 

l2S] 

1/2 

3s23p4  ( iP)4s 

l*P] 

5/2,  3/2,  1/2 

l2P] 

3/2,  1/2 

( lD)4s' 

[\D\ 

3/2,  5/2 

OS)4s’ 

[2S] 

1/2 

3s23p4  ( 3P)3d 

[*D] 

5/2,  3/2,  1/2 

[*/] 

5/2,  3/2 

l2P] 

1/2,  3/2 

[V] 

1/2,  3/2,  5/2 

l2P] 

5/2 

( lD)3d' 

t201 

3/2,  5/2 

l2P] 

5/2 
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>  Odd  parity  energy  levels  with  excited  electron  configurations 
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electron  impact  excitation: 

radiative  decay: 

electron  impact  ionization: 
radiative  recombination: 

three  body  recombination: 


Arl  and  Aril  can  be  excited  from  one  energy  level  to  a 
higher  energy  level  due  to  inelastic  scattering  collisions  with 
incident  electrons. 

Excited  levels  of  Arl  and  Aril  can  spontaneously  decay  to 
lower  energy  levels. 

Arl  can  be  ionized  by  electron  impact  creating  Aril. 

Aril  can  recombine  with  secondary  electrons  creating  Arl 
and  a  photon. 

Aril  can  recombine  with  secondary  electrons  in  the 
presence  of  Arl. 


The  following  assumptions  were  made  to  restrict  the  model: 


a)  All  recombinations  result  in  ground  level  neutral  argon. 

b)  Arl  that  experience  ionization  produce  Aril  in  either  the  first  two  energy  levels,  i.e., 
those  with  the  (Ne)  3s23ps  electron  configuration  or  the  excited  energy  levels  with 
the  (Ne)  3s23p44p  electron  configuration. 

c)  Aril  ions  that  are  produced  in  the  two  energy  levels  of  the  ground  electron  configu¬ 
ration  are  generated  with  respective  probabilities  PI  and  P2  such  that  PI  +  P2  =  1. 
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Since  these  two  levels  are  so  closely  spaced,  the  distribution  of  ions  created  in  each 
of  these  two  energy  levels  were  assumed  to  be  in  proportiooio  their  degeneracies,  i.e.. 
Pi  =  J  and  P2  =  $. 

d)  Since  the  energy  levels  corresponding  to  the  electron  configuration  (Ne)  3s23p*4p 
are  the  highest  energy  configuration  under  consideration,  these  4 p  levels  will  not 
experience  excitation  to  higher  levels  nor  will  they  have  gain  terms  due  to  the  higher 
levels  decaying  down  to  this  configuration. 

e)  There  is  a  small  number,  on  the  order  of  one  percent,  of  neutral  Arl  atoms  which 
ionize  directly  to  Aril*  in  excited  energy  levels  with  the  (Ni)  3s23p*4p  configu¬ 
ration. 


3.  Population  rate  equations 

The  rate  equations  for  the  population  density  of  the  /th  energy  level  can  be  expressed  in 
the  general  form: 


dN(D 

dt 


=  G  —  L 


(3) 


where  G  is  the  gain  term  and  L  is  the  loss  term  for  N{1),  the  population  density  of  energy 
level  /.  Given  the  model  and  assumptions  discussed  in  the  previous  section,  the  rate  equa¬ 
tions  were  written  for  the  reactions  shown  in  table  2. 

The  rate  equation  for  neutral  argon  with  the  ground  level  electron  configuration  is: 

=  2  A,(*./)N,(M)  -  f  2  2 

"*  *>l  1*1-1  M-fl  +  l  J 

-  2 *, (/,*)*!</,/)  +  2  U<*>  +  Ms)  2*«i<'.')]n,(*,o  (4) 

><i  *-i  L  i  J 

and  the  rate  equation  for  levels  of  neutral  argon  with  excited  electron  configurations  is: 


Tabu  2.  Chemical  reactions  included  in  the  model 


1)  Electron  excitation  of  Arl: 

2)  Electron  ionization  of  Arl: 

3)  Radiative  decay  of  Arl: 

4)  Direct  excitation  of  Arl: 

3)  Electron  excitation  of  Aril: 

6)  Radiative  decay  of  Aril: 

7)  Radiative  recombination: 

8)  3  body  recombination: 

9)  Creation  of  Arl(meta): 

10)  Ionization  of  Arl(meta): 


e  +  Arl  -»  Arl*  +  e 
e  +  Arl  -*  ArII(3p)  +  2e 
Arl*  -» Arl  +  hr 
e  +  Arl  -*  Aril*  (<p)  +  2e 
e  +  Aril  -*  Aril*  +  e 
e  +  Aril*  -  Aril*  +  e 
ArII*(4p)  -•  Aril*  +  hr 
Aril*  -» Aril  +  hr 
e  +  Aril  -» Arl  +  hr 
e  +  Aril*  Arl  +  hr 
e  +  Aril  +  Arl  -*  Arl  +  Arl 
e  +  Arl  -•  Arfmeta)  +  e 
e  +  Arl(meu)  -•  Aril  +  Xe 


Jt 
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Sll 


dNiU) 

dt 


J}ktVJ)N,{i.t) 

i<  1 


^  2  A, </,/)  + 


E 


*ic  \ 


»t 

E  <*il>8u82)NAgut) 

*2=*!+' 


N,(/./) 


(5) 


k\(i,k) 

Ai(k,i) 

aU.gi.gi) 

0r(g) 

03(8) 

K(g,r) 

N,(/,0 


is  the  pumping  rate  {sec-1]  for  generating  excitations  from  energy  level  /  to 
energy  level  k  in  Arl  due  to  electron  impact  excitation  by  the  beam  electrons 
and  all  secondary  electrons  with  energy  above  the  threshold  £*; 
is  the  Einstein  A  coefficient  [sec-1]  for  radiative  decay  from  Arl  in  energy 
level  k  to  Arl  in  energy  level  <  due  to  the  spontaneous  emission  of  radiation; 
is  the  ionization  coefficient  [cmVsec]  for  electrons  in  group  gi  causing  an 
ionization  of  Arl  in  energy  level  /  resulting  in  an  emitted  electron  in  energy 
group  g2; 

is  the  radiative  recombination  coefficient  [cmVsecJ  for  electrons  in  group 
g  recombining  with  Aril  to  create  Arl; 

is  the  three  body  recombination  coefficient  IcmVsecl  for  electrons  in 
group  g  recombining  with  Aril  to  create  Arl; 

is  the  number  density  [cm-3]  of  electrons  with  kinetic  energies  in  the 
interval  (Tt<T<  7"f_i)  at  time  t;  and 

is  the  number  density  [cm-3]  of  Arl  atoms  in  energy  level  /  at  time  t. 


The  rate  equation  for  the  two  ground  levels  of  singly  ionized  argon  is: 


dN,U\  22  N*  N‘ 

at  /-i  *i-i  *2-*(+i 

-  S/MfWU.OAM/.O 

k> I  *-l 

Nt 

+  2  Au(k,l)Nu(k,t)  -  2  03 (g)Ne(g.n  (6) 

*>1  *«1  !' 

The  rate  equation  for  levels  of  singly  ionized  argon  with  (Ne)  3s23p44sor  (Ne)  3s23p43d 
electron  configurations  is: 


N  r  22 

^4 —  =  H/cu(U)N»(itt)  -  s  U<*)  +  03(g)  2  Nx(l\t)  Ne(g,t)N„(l,t) 
at  ,<l  *.i  i  r-i 

-  2  *»</.*)  +  2^u(/.o  |a^i.(/.o  +  2'4u<*.wv„(*,o.  (7) 

.*>1  «l  J  *>l 

The  rate  equation  for  levels  of  singly  ionized  argon  with  (Ne)  3s23p44p  electron  config¬ 

uration  is 


=  2*n(/./)Arn(/.0  -  ^  U(g)  +  03(g)  2  W,t)\N.(g,t)*n0.t) 
at  /<i  *-i  l  r-i  J 

-'ZA„(U)N„(l,t)+  2  E  E  /(/)«(/'.*i.fc)N*(gi.ONi(/',0 

/<!  rm |  *,-l  *2-*!+l 

(8) 
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where 

Arn (/.At)  is  the  pumping  rate  [sec-1]  for  generating  excitations  from  energy  level  r 
to  energy  level  k  in  Aril  due  to  electron  impact  excitation  by  the  beam 
electrons  and  all  secondary  electrons  with  energy  above  the  threshold  £*; 

*ihtrm  (»’.*)  is  the  pumping  rate  [sec-1]  for  generating  excitations  from  energy  level  / 
to  energy  level  k  in  Aril  due  to  electron  impact  excitation  by  the  thermal 
electrons  with  energy  above  the  threshold  £»; 

A, ,(£,/)  is  the  Einstein  A  coefficient  [sec-1]  for  radiative  decay  from  Aril  in  en¬ 
ergy  level  k  to  Aril  in  energy  level  i  due  to  the  spontaneous  emission  of 
radiation; 

Nn(l,t)  is  the  number  density  [cm-3]  of  Aril  ions  in  energy  level  /  at  time  /;  and 

/(/)  is  the  fraction  of  Arl  directly  ionized  to  the  allowed  3s23p*4p  energy 

levels. 

The  pumping  rates  kx,ku  [sec-1]  for  exciting  Arl  and  Aril  from  energy  level  i  to  energy 
level  k  are  of  the  following  form  (Bretagne  et  al.  1988): 


kik(T„)  = 


Vik(Tb)Lj{Tb) 


where  L,(Tb)  is  the  total  loss  function  for  an  electron  with  kinetic  energy  Tb  incident  upon 
a  target  in  energy  level  r,  Jb  is  the  current  density  of  the  electron  beam,  e  is  the  charge  of 
an  electron,  itut(Tb)  is  the  efficiency  for  an  electron  with  kinetic  energy  Tb  generating  a 
transition  in  the  target  from  energy  level  r  to  energy  level  k,  and  Ea  is  the  energy’  differ¬ 
ence  between  levels  k  and  /. 

The  ionization  coefficients  [sec-1]  are  calculated  similarly: 


k0(Tb)  = 


nu(Tb)L,(Tb) 


where  ij,*(r6)  is  the  efficiency  for  an  electron  with  kinetic  energy  Tb  generating  an  ioniza¬ 
tion  in  the yth  shell  of  a  target  in  energy  level  /,  and  /</  is  the  ionization  potential  of  the  yth 
shell  of  a  target  in  energy  level  i.  Lj(Tb),  Jb,  and  e  are  defined  the  same  as  for  (9). 

The  cascade  electrons  experience  elastic  electron-electron  collisions  which  result  in  the 
thermalization  of  the  electrons  within  each  energy  group.  The  time  for  this  process  to  take 
place  is  on  the  order  of  a  nanosecond.  Since  some  electrons  in  energy  bins  below  the  in¬ 
elastic  threshold  will  be  scattered  elastically  to  energies  above  the  threshold,  their  contri¬ 
bution  must  also  be  considered  even  though  the  actual  elastic  scattering  will  not  be  performed 
in  the  simulation.  Boltzmann  plots  of  selected  experimentally  measured  intensities  of  spec¬ 
tral  emissions,  believed  to  be  dominated  by  the  plasma  electrons,  indicated  an  electron  tem¬ 
perature  of  roughly  1.7  eV  (Brake  &  Repetti  1988).  This  is  also  the  temperature  which 
characterizes  the  Maxwellian  distribution  that  approximates  the  predicted  distribution  of 
electrons  after  elastic  scattering.  For  each  transition  from  energy  level  i  to  energy  level  k, 
an  excitation  rate  will  be  calculated  by  evaluating  the  following  at  T  »  1.7  eV: 


*«h«rm»i(».*)  =n(a,kv)  =  f  Fmtx(.kBT)otk(E  =  kbT)v(E  =  kbT)dE.  (11) 

This  thermal  rate  Ar(hmMj(i, k)  in  [sec-1]  was  then  added  to  the  term  k(i,k)  previously  cal¬ 
culated  for  the  rate  of  transitions  created  by  the  beam  electron  and  secondary  cascade  elec¬ 
trons  alone  to  give  a  total  excitation  rate. 

It  was  also  necessary  to  write  rate  equations  for  energy  groups  of  secondary  electrons 
since  the  recombination  coefficients  are  functions  of  electron  energy.  The  primary  electrons 
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represent  a  special  case  of  electrons  all  with  kinetic  energy  equal  to  the  beam  electron  en¬ 
ergy,  300  keV  in  this  case.  The  beam  electron  rate  equation  is  deduced  from  the  measured 
current,  the  dimensions  of  the  cathode,  and  the  energy  of  the  beam  where 


dN^jl,  t)  = 

dt 


C, 


t<  r, 


0  T,  <  /  <  T2 
~C2  T2  £  /  <  Tj 
0  t>r3 


(12) 


with  Ci  =  8.2  x  10,7(cm-3  s-1].  C2  =  2.7  x  10,7[cm"3  s-,J.  ri  =  25  ns,  t2  =  225  ns,  and 
rj  =  300  ns  to  match  the  experimental  conditions  (Brake  et  al.  1986;  Brake  &  Repetti 
1988). 

The  energy  bins  for  secondary  electrons  encompassed  electrons  with  energies  of  200, 100, 
50, 20,  10,  5, 2,  1, 0.5, 0.2,  and  0.1  eV.  The  upper  limit  was  chosen  from  the  examination 
of  the  calculated  distribution  of  secondary  electrons  after  one  ionization  of  a  valence  elec¬ 
tron  from  ground  level  neutral  argon.  The  number  of  secondary  electrons  with  kinetic  ener¬ 
gies  over  200  eV  was  negligible.  An  energy  bin  “g"  is  defined  such  that  it  contains  all 
electrons  with  kinetic  energies  Tt<T<  Tt.,.  The  gain  term  for  any  given  energy  bin  in 
the  secondary  electron  rate  equation  is  due  to  electrons  of  higher  energy  (created  from  ioniz¬ 
ing  collisions)  that  scatter  down  into  the  energy  bin.  The  loss  term  is  due  to  electrons  which 
suffer  scattering  collisions  which  cause  them  to  lose  sufficient  energy  to  drop  down  to  a 
lower  energy  bin  as  well  as  electrons  in  the  bin  which  experience  recombination  with  an  Aril 
ion. 

The  rate  equation  for  secondary  electrons  of  group  g  is  given  by  the  following  expression: 

45  a(/.G«.*)N,(/.ON,(G„0 
dt  f.i  c,-i 

+  §  ?>xU,C2,g)NtU,t)NAGt,t) 

Oi-l  /-i 

-  L  2x</.2,GI)N,(/.r)N,(*.0 

Oft  /-l 

[22  1  » 

0,(g)  +  ft(g)  s  N, </',*)  K(8.0  LNii</.0 
r-i  J  /-i 

where  a  and  P  are  the  ionization  and  recombination  coefficients,  as  previously  defined.  The 
rate  [cm3  sec-1]  of  creation  of  electrons  due  to  ionization  is  calculated  as  follows: 

<x(Tt)  =  o^lTJv,  (14) 

which  can  be  rewritten  as 

f  <foioa(7f>£)  Jr- 
°<rt)  =  vtJ  -  dE 


(13) 


(15) 


by  making  the  following  approximation 


f  ffejg  ^TvEidE> 
J  dE 


j* Edoiea(TrE)  d£ 

J _ _ 

jdE 


3 7- 


(16) 
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and  using  the  definition  of  the  loss  function  given  in  McGarrah  &  Brake  (1990),  the  ion¬ 
ization  coefficient  can  be  expressed  as 


AE 


(17) 


The  scattering  coefficients  x(l.gt>gi)  (cm1  sec-1]  represent  the  rate  at  which  electrons 
experience  scattering  collisions,  namely,  excitation  or  ionization,  that  result  in  electrons  hav¬ 
ing  energy  characteristics  of  one  energy  group  to  an  energy  within  the  interval  of  a  lower 
energy  group.  The  scattering  (excitation  and  ionization)  coefficient  can  be  determined  using 
the  total,  as  opposed  to  ionization,  loss  function. 


X(7i)  = 


^Tor(Tt)vt 

AE 


(18) 


where  for  energy  group  g,  in  both  cases, 


AE=Tt.,  -  Tt.  (19) 

The  series  of  coupled  rate  equations  was  solved  using  a  software  package  of  FORTRAN 
codes  created  at  Lawrence  Livermore  Laboratory  known  as  EPISODE  (experimental  pack¬ 
age  for  integration  of  systems  of  ordinary  differential  equations,  see  Hindmarsh  and  Byrne 
(1975)). 


4.  Results 

The  electron  density  can  be  obtained  by  summing  each  of  the  calculated  secondary  elec¬ 
tron  energy  groups  N,(g,t)  over  all  g.  The  result  of  these  calculations  at  a  pressure  of  50 
Torr  is  shown  in  figure  1.  The  plasma  density  is  roughly  linear  over  the  period  of  the  beam 
pulse  and  then  remains  constant.  This  behavior  is  in  agreement  with  literature  (Miller  1982) 
which  indicates  that  the  partial  time  derivative  of  the  plasma  density  due  to  the  beam  dec- 
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Figure  1.  Plasma  density  versus  time  at  pressure  equal  to  SO  Torr. 
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irons  is  proportional  to  the  beam  electron  density.  In  this  case  the  beam  a 
for  20 0  nanoseconds,  suggesting  a  linear  time  dependence  for  the  plasm  .cnst(*  *s  con<i,anf 
will  eventually  decrease  due  to  recombination  long  after  the  beam  pulse  s-SU>  va*uc 
was  calculated  for  pressures  up  to  750  Torr  with  peak  density  of  1.7  x  ionllar  _^hav'or 
Figure  2  shows  the  time  dependent  populations  of  9  upper  levels  of  the  10 ■ 
the  visible  range  as  predicted  by  the  simulation.  They  are  labeled  4p(i)  through^ " t qTk!.0 
cause  they  all  have  (Ne)  3s23p44p  electron  configuration.  The  remainder  of  their  des'^  ** 
tion  is  indicated  in  table  3.  lgna‘ 

Three  of  the  line  intensities  were  dramatically  larger  than  the  other  seven.  These  three 
emissions  have  upper  levels  with  4 p'  and  4 p"  valence  electrons  (levels  7,  8,  and  9  of  ta¬ 
ble  3),  i.e.t  they  are  among  the  highest  energy  levels  observed.  A  thermal  model  would 
indicate  that  these  populations  should  be  small  relative  to  the  less  energetic  levels  with  4 p 
valence  electrons.  The  mechanism  of  direct  ionization  indicated  by  Reaction  (4)  in  table  2 
must  satisfy  the  selection  rules  (Koozekanani  1966)  of  equations  (1)  and  (2);  therefore  it  can 
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Figure  2.  Populations  of  the  nine  (Ne)  3s23p44p  excited  electron  energy  levels  of  Aril  under  con¬ 
sideration  at  pressure  equal  to  100  Torr.  See  table  3  for  level  configuration. 


Table  3.  Identification  of  energy  levels  for  figure  2 

Label  Configuration-Term  J 

4p(l)  3s*3p*  (*P)4p  [*D]  1/2 

4p<2)  (JP)4p  [2P]  1/2 

4p(3)  ( 5P)4p  l2PJ  3/2 

4p(4)  (3P)4p  (4P1  3/2 

4p(5)  ixD)4p'\\D)  5/2 

4p(6)  (lD)4p'[lF]  7/2 

4p(7)  ( ,D)4p'  1 2PJ  3/2 

4p(8)  ( *S)  4p*  [  2PJ  1/2 

4p(9)  ( lS)  4p'[2P]  3/2 
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only  contribute  to  the  populations  of  the  low  energy  4p(l)  through  4p(4)  levels  and  the 
4p'  level  with  J  =  $  (i.e.,  4/>(7))  and  the  two  4 p"  levels  (4/>(8)  and  4p(9)).  Due  to  this 
direct  ionization,  the  two  4 p"  energy  levels  and  the  4 p'  [2/*3/2l  energy  level  are  substan¬ 
tially  larger  than  the  other  energy  levels  and  appear  to  rise  sharply  with  the  incidence  of 
the  electron  beam.  They  remain  flat  while  the  beam  pulse  is  constant  and  rapidly  decline 
as  it  is  turned  off.  The  smallest  populations,  corresponding  to  lower  energy  levels,  increase 
slowly  with  time  or  remain  roughly  constant  even  after  the  beam  pulse.  This  is  similar  to 
the  time  dependence  of  the  plasma  electron  density. 

The  time-dependent  populations  of  the  upper  levels  were  used  to  generate  relative  inten¬ 
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Figure  3a.  Comparison  of  normalized  calculated  intensities  (solid  boxes)  with  measured  values  (open 
boxes)  for  10  experimentally  observed  wavelengths  at  100  Torr. 
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Figure  3b.  Comparison  of  normalized  calculated  intensities  (solid  boxes)  with  measured  values  (open 
boxes)  for  10  experimentally  observed  wavelengths  at  200  Torr. 
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sity  values  so  that  they  could  be  compared  to  experimental  results  obtained  from  an  intense, 
relativistic  electron  beam  (300  keV,  lkA,  300  ns)  injected  into  a  vacuum  vessel  which  was 
back  filled  with  argon  at  pressures  ranging  from  1  to  750  Torr  (Brake  et  at.  1986;  Brake 
&  Repetti  1988).  Since  relative  rather  than  absolute  intensities  were  measured  in  the  exper¬ 
iment,  the  model  intensities  were  scaled  to  the  experimentally  observed  intensities  by  nor¬ 
malizing  them  to  the  4053  A  line. 

The  results  for  pressures  from  100  Torr  to  750  Torr  are  given  in  figure  3a-3d  and  show 
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Figure  3c.  Comparison  of  normalized  calculated  intensities  (solid  boxes)  with  measured  values  (open 
boxes)  for  10  experimentally  observed  wavelengths  at  380  Torr. 
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figure  3d.  Comparison  of  normalized  calculated  intensities  (solid  boxes)  with  measured  values  (open 
boxes)  for  10  experimentally  observed  wavelengths  at  750  Torr. 
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excellent  agreement.  In  figure  4a  and  4b  the  calculated  ratio  of  the  intensity  of  the  4053  A 
line  to  the  two  other  enhanced  lines  at  3995  A  and  4278  A  were  compared  to  the  experi¬ 
mentally  observed  values  over  a  range  of  pressures.  The  calculated  values  agree  extremely 
well  with  measured  values  and  also  support  the  concept  that  these  lines  are  dominated  by 
beam  electron  densities,  constant  over  pressure,  and  Aril  densities,  linear  with  pressure. 
If  all  three  intensities  are  linear  in  pressure,  their  ratios  will  then  be  constant  with  pressure. 
The  intensities  of  the  three  lines  at  3995,  4053,  and  4278  A  are  presented  as  a  function  of 
pressure  in  figure  5.  The  calculated  values  are  indeed  linear  with  pressure.  The  experimental 
values  concur  at  low  pressures  but  begin  to  diminish  at  higher  pressures  where  the  analyt¬ 
ical  assumption  of  an  optically  thin  plasma  may  begin  to  break  down. 

The  high  energy  4 p"  and  4p'[2P3/2]  energy  levels  can  be  reached  due  to  direct  ioniza¬ 
tion  of  neutral  argon  in  the  ground  level.  This  is  possible  because  the  transition  from  Arl 
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Figure  4a.  Comparison  of  the  calculated  (solid  line)  and  measured  ratios  of  the  normalized  inten¬ 
sity  for  the  4053  A  to  the  3995  A  line. 
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Figure  4b.  Comparison  of  the  calculated  (solid  line)  and  measured  ratios  of  the  normalized  inten¬ 
sity  for  the  4053  A  line  to  the  4278  A  line. 
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Figure  5.  Comparison  of  calculated  (solid  line)  and  experimental  intensities  of  3995  A  (■),  4053  A 
(•).  and  4278  A  (a)  lines  versus  pressure. 
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to  those  levels  still  satisfies  equation  (1).  Although  only  about  19b  of  ionization  is  believed 
to  produce  excited  energy  levels  (Bennett  et  al.  1966),  this  still  provides  a  mechanism  to 
populate  these  high  energy  levels  which  will  dominate  over  the  three  step  procedure: 

e  +  Arl- 


ArII(3p)  +  2e 
e  +  ArII(3p)  -» Aril  (4s,  3d)  +  e 
e  +  ArII(4s,3d)  -►  Arll(4p)  +  e. 


(20a) 

(20b) 

(20c) 


The  beam  and  cascade  electrons  are  able  to  drive  this  process  because  the  direct  ioniza¬ 
tion  of  Arl  to  energetic  Aril*  represents  a  transition  of  over  30  eV  and  is  therefore  unlikely 
to  be  generated  by  the  thermal  electrons  characterized  by  an  energy  of  less  than  2  eV.  How¬ 
ever,  direct  ionization  would  be  a  more  probable  transition  for  the  cascade  electrons  with 
energies  on  the  order  of  100  eV.  The  significance  of  this  effect  is  indicated  by  figure  6,  where 
the  results  of  a  simulation  performed  without  the  inclusion  of  direct  ionization  are  com¬ 
pared  with  experimental  intensities  for  the  10  observed  transitions,  again  scaled  to  the  4053 
A  line.  The  results  are  seen  to  differ  dramatically. 
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5.  Conclusions 

The  solution  of  coupled  differentia]  rate  equations  describing  the  chemical  kinetics  of  an 
intense  relativistic  electron  beam  interacting  with  neutral  argon  indicates  that  the  gas  ex¬ 
periences  ionization  of  about  19b  of  the  argon.  The  plasma  density  increases  linearly 
throughout  the  beam  pulse  and  remains  for  some  time  after  the  beam  is  turned  off.  Low 
energy  excited  levels  of  Aril  have  time  profiles  similar  to  the  plasma  density  and  appear 
to  be  generated  by  the  thermal  electrons.  High  energy  4p'  levels  and  the  4p'  level  that  can 
be  reached  by  direct  ionization  of  Arl  have  the  same  temporal  profile  as  the  electron  beam 
current.  These  levels  are  populated  almost  entirely  by  direct  ionization  by  the  beam  and  cas¬ 
cade  electrons,  allowing  them  to  be  more  densely  populated  than  the  lower  energy  levels. 
The  high  energy  Ap"  levels  then  have  a  large  Einstein  A  coefficient  for  decaying  to  4s'  en- 
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Wavelength  (A) 

Figure  6.  Predicted  relative  intensities  (dashed  box)  for  a  simulation  which  does  not  include  direct 
excitation  from  Arl  to  4p  levels  of  Aril  compared  with  the  experimental  values  (solid  box). 


ergy  levels,  which  in  turn  have  high  beam  deposition  efficiencies  for  being  pumped  by  the 
electron  beam  back  up  to  the  4 p"  levels.  In  this  way  certain  transitions  between  energetic 
excited  levels  of  Aril  are  pumped  by  the  electron  beam.  These  transitions  between  high  en¬ 
ergy  excited  levels  would  be  negligible  if  the  mechanism  of  direct  ionization  was  not  con¬ 
sidered.  The  model  achieved  excellent  agreement  with  the  experimental  measurements  of 
relative  intensities  of  spectral  emissions  from  argon  ions  over  a  range  of  pressures  from  SO 
to  750  Torr  and  should  prove  useful  for  other  studies  of  beam  pumped  rare  gases. 
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